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Abstract 

In the recent publication (Journal of Geometry and Physics, 33 (2000) 
23-102 ) we have demonstrated that dynamics of 2+1 gravity can be de- 
scribed in terms of train tracks. Train tracks were introduced by Thurston 
in connection with description of dynamics of surface automorphisms. In 
this work we provide an example of utilization of general formalism devel- 
oped earlier. The complete exact solution of the model problem describing 
equilibrium dynamics of train tracks on the punctured torus is obtained. 
Being guided by similarities between the dynamics of 2d liquid crystals 
and 2+1 gravity the partition function for gravity is mapped into that 
for the Farey spin chain. The Farey spin chain partition function, fortu- 
nately, is known exactly and has been thoroughly investigated recently. 
Accordingly, the transition between the pseudo-Anosov and the periodic 
dynamic regime (in Thurston's terminology) in the case of gravity is being 
reinterpreted in terms of phase transitions in the Farey spin chain whose 
partition function is just the ratio of two Riemann zeta functions. The 
mapping into the spin chain is facilitated by recognition of a special role 
of the Alexander polynomial for knots/links in study of dynamics of self 
homeomorphisms of surfaces. At the end of paper, using some facts from 
the theory of arithmetic hyperbolic 3-manifolds ( initiated by Bianchi in 
1892), we develop systematic extension of the obtained results to noncom- 
pact Riemann surfaces of higher genus. Some of the obtained results are 
also useful for 3+1 gravity. In particular, using the theorem of Margulis, 
we provide new reasons for the black hole existence in the Universe: black 
holes make our Universe arithmetic. That is the discrete Lie groups of 
motion are arithmetic. 
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1 Introduction and summary 



The Riemann zeta function ((/3) has been an object of intensive study in 
both mathematics [1-3] and physics [4,5] for quite some time. The reason for 
physicists interest in this function can be easily understood if one writes it in 
the form of a partition function Z(/3) given by 



If (3 is interpreted as the inverse temperature then, naturally, questions arise: 

a) what is the explicit form of the quantum mechanical Hamiltonian whose 
eigenvalues E n are given by E n = In n ? 

b) can such system undergo phase transition(s) if one varies the temperature? 
The goal of providing answers to both questions is at the forefront of current 
research activities both in physics [4,5] and mathematics [6]. Answers to these 
questions are being sought in connection with theories of random matrices and 
quantum chaos [4,5], non-commutative geometry [6] and Yang-Lee zeros [7]. 
According to the theory of Yang and and Lee the problem of existence of phase 
transitions can be reduced to the problem of existence of zeros of the partition 
function in the complex z-planc (where z may be related to either fugacity or 
the magnetic field, etc.). In the case of Z((3), Eq.(l.l), one is also looking at 
analytic behavior of the partition function in the complex j3— plane. Riemann 
had conjectured that, 



provided that Re t m ^ for all integer m's. This conjecture is known in the 
literature as Riemann hypothesis. Stated differently, the Riemann hypothesis is 
equivalent to the statement that all " nontrivial" zeros of the partition function 



for m=l,2,... and B2 m being the Bernoulli numbers. 

Combinations of the Riemann zeta functions are also of physical interest. In 
particular, in this paper we shall be concerned with the following combination 



oo 




(1.1) 




(1.2) 





(1.4) 



where <p(n) is the Eulcr totient function, 



<j)(n) 




(1.5) 
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which is just the number of numbers less than n and prime to n, provided that 
n = p™ 1 ■ ■ ■ p™ T , and pi, etc. are primes with respect to n. This partition 
function had appeared in mathematical physics literature in connection with 
the partition function for the number-theoretic spin chain [8-10] and in connec- 
tion with calculations of the scattering S-matrix for the " leaky torus" quantum 
mechanical problem [11-13]. Remarkably enough, the results for the number- 
theoretic spin chain can be obtained as well from earlier works on mode locking 
and circle maps [14,15], e.g. see in particular Eq.(30) of Ref.[15], as the authors 
of Ref.[10] acknowledge. This fact is not totally coincidental as we shall explain 
below (in sections 2 and 3) and has been already anticipated based on our earlier 
works [16,17] on dynamics of 2+1 gravity. 

In this paper we would like to demonstrate that the partition function, 
Eq.(1.4), can adequately describe statistical mechanics of Einstcinian 2+1 grav- 
ity if the underlying surface is the punctured torus. The restriction to the 
punctured torus case is not too severe and is motivated mainly by illustrative 
purposes: recall, that both the Seifert surfaces of the figure eight and the trefoil 
knots are just punctured toruses [18]. This observation allows us to make an 
easy connection between the dynamics of surface self homeomorphisms and the 
associated with its time evolution 3-manifolds which fiber over the circle [17] . 
These manifolds are just complements of the figure eight and trefoil knots in 
S 3 respectively.The Seifert surfaces of more complicated knots may naturally 
be of higher genus but, since both the trefoil and the figure knots belong to 
the category of fibered knots, only those knots and links which are fibered and 
the associated with them Seifert surfaces are relevant to the dynamics of 2+1 
gravity [17]. The 3- manifolds associated with the figure eight and the trefoil 
knots are fundamentally different: the first one is known to belong to the sim- 
plest representative of the hyperbolic manifolds while the second corresponds to 
the so called Seifert fibered manifolds [19]. The surface dynamics associated with 
the first is associated with pseudo-Anosov type of surface self homeomorphisms 
while the second one is associated with periodic self homeomorphisms. Both 
types of 3-manifolds are topologically very interesting and potentially contain 
wealth of useful physical information. In this work we only initiate their study 
with hope of returning to this subject in future publications. 

In order for the reader to keep focus primarily on physical aspects of the 
problem, we feel, that some simple explanation of what follows is appropri- 
ate at this point. To avoid repetitions, we expect that our readers have some 
background knowledge of the results presented in our earlier published papers 
[16,17]. In particular, to help our intuition, we would like to exploit the fact 
that statics and dynamics of 2+1 gravity is isomorphic with statics and dy- 
namics of textures in two dimensional liquid crystals. According to the existing 
literature on liquid crystals, e.g. [20], the liquid crystalline state can be found 
in several phases which physicists classify as liquid, solid, hexatic and gas. In 
mathematical literature the textures, e.g. like those in liquid crystals, are known 
as foliations [19,21]. Dynamics of textures is known accordingly as dynamics 
of foliations. Some of these foliations may contain singularities. These singular- 
ities are mistakenly being treated as Coulombic charges (while in 2+1 gravity 
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is is well documented [17,22] that these singularities do not interact) since the 
nonorientable line fields arc being confused with the orientable vector fields. 
The phase transitions in two dimensional liquid crystals are described in terms 
of the phase transitions in 2 dimensional Coulomb gas. These are known as 
the Kosterlitz-Thouless type transitions [23]. Although mathematically such 
explanation of phase transition is not satisfactory, nevertheless, one has to re- 
spect the experimental data associated with the liquid crystalline phases. The 
extent to which the Kosterlitz-Thouless intcrprcpretation of phase transition is 
appropriate is discussed in the Appendix A.l. where it is being argued that, 
by analogy with transitions in the liquid helium (where the Kosterlitz-Thouless 
interpretation is normally used) the dynamical phase transition in 2+1 gravity 
resemble to some extent the Bose gas condensation type of transition. This anal- 
ogy is incomplete, however, and is only being used for the sake of comparison 
with the existing literature. As Remark 4.1. indicates, the partition function 
of 2+1 gravity without any appriximations can be recast into the Lee- Yang 
form [7]. The calculations associated with such form require knowledge of dis- 
tribution of zeros of the Riemann Zeta function and, hence, effectively, the proof 
of the Riemann hypothesis. Since this proof is not yet available (see, however, 
Ref.[5]) we employ alternative methods associated with recently developed ther- 
modynamic/statistical mechanic formalism for description of phase transitions 
in the number-theoretic Farey spin chains [8-10]. To prove that dynamical tran- 
sitions in gravity can be described in terms of transitions in the Farey spin 
chains several steps are required. In our earlier works [16,17] we have demon- 
strated that dynamic of 2+1 gravity is best described in terms of dynamic of 
train tracks. In section 2 we demonstrate how dynamic of train tracks can be 
mapped into dynamic of geodesic laminations. In turn, the dynamic of geodesic 
laminations is reformulated in terms of the sequence of Dchn twists. This se- 
quence is actually responsible for the dynamic in the Teichmuller space 
of the punctured torus. Such dynamic is subject to the number-theoretic con- 
straints associated with the Markov triples. The Markov triples had been known 
in physics literature for a while in connection with the trace maps [24] used 
for description of quasicrystals, Id tight band Scrodinger equations, etc. In 
this work the Markov triples play a somewhat different role: they make the 
set of closed nonpcriphcral geodesies on the punctured torus discrete. Math- 
ematically, the sequence of Dehn twists is written in terms of the product of 
the "right" and the "left" 2x2 Dchn matrices. The modulus of eigenvalues(the 
stretch factors) associated with such matrix product can be either greater than 
one or equal to one. In the first case one is dealing with the pseudo-Anosov 
and in the second, with the periodic (Seifert-fibered) dynamical regime. The 
results of section 2 (for the figure 8 and the trefoil knots) acquire new meaning 
in section 3 where they are reobtained with help of the associated Alexander 
polynomials. The stretch factors of section 2 are reobtained as zeros of the re- 
lated Alexander polynomials. In the same section 3 we discuss the fiber bundle 
construction of 3-manifolds complementary to the figure eight and the trefoil 
knots in S 3 . The sequence of Dehn twists, discussed in section 2, in this section 
is being associated with the operation of Dehn surgery (Dchn filing) performed 
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on the 3-manifold related to the figure eight knot. The stretch factors produced 
as result of such surgery are reobtained with help of the Mahler measures which 
allow us to reinterpret these factors in terms of the topological entropies. Intro- 
duction of the Mahler measures, in addition, allows us to make direct connection 
between the dynamical phase transitions and the thermodynamic transitions in 
the sense of Yang and Lee [7] : zeros of the Alexander polynomial play similar 
role in dynamics as Yang-Lee zeros in thermodynamics. In section 4 we pro- 
vide direct connection between the results of section 3 just described and the 
statistical mechanics formalism developed for the number-theoretic Farey spin 
chains [8-10]. With such connection established, dynamical phase transitions in 
2 +1 gravity can be treated with formalism which is more familiar to physicists. 
Unfortunately, this more familiar formalism is applicable (at least at the present 
level of our understanding) only to the case of punctured torus. Fortunately, 
the final results obtained with help of such type of formalism can be reobtained 
in several different and independent ways. We discuss these alternative ways at 
the end of section 4 and in section 5 which is entirely devoted to development 
and refinements of these alternatives. The major reason of doing this lies in 
the opportunity of extension of the punctured torus results to the noncompact 
Ricmann surfaces of any genus. This is accomplished using some results from 
the scattering theory for Poincare' (Eisenstein) series acting on 3-manifolds. 
These series had been recently discussed in our earlier work, Rcf.[25], to which 
we refer for more details. The key result of section 4, the partition function for 
the case of punctured torus, happens to coincide with the scattering S-matrix 
(up to unimportant constant) obtained for 3-manifolds with one Z®Z cusp [26]. 
3-manifolds containing multiple Z+Z cusps are associated with fibered links (as 
explained in section 3 and the Appendix A. 3). The S-matrix for this case has 
been also obtained recently [27]. The determinant of this matrix produces the 
desired exact partition function for 2+1 gravity. In addition, the formalism al- 
lows us to obtain the volumes of the associated 3-manifolds exactly. Extension 
of the punctured torus results to surfaces of higher genus requires some careful 
analysis of the nontrivial mathematical problem of circle packing. This fact has 
some profound impact at development of the whole formalism. It happens, that 
such type of problems had been comprehensively studied by Bianchi already in 
1892 [28] who championed study of the arithmetic hyperbolic manifolds. The 
notion of arithmcticity is rather involved. Since in physics literature (to our 
knowlcgc) it did not find its place yet, in Appendices A. 2 and A. 3 we supply 
the essentials needed for the uninterrupted reading of the main text. Surely, 
selection of the material in these appendices is subjective. But, it is hoped, that 
interested reader will be able to restore the missing details if it is required. The 
arithmcticity of 3-manifolds associated with 2+1 gravity stems from some very 
deep results of Riley [29] and Margulis [30] which we discuss to some extent in 
Appendix A. 3. The arithmeticity leads to some restrictions on groups of mo- 
tions in symmetric spaces (e.g. hyperbolic space is symmetric space). Thanks to 
the Margulis Theorem A. 3. 11. and some results of Hclgason [31] and Besse [32], 
the notion of arithmeticity is extendable to 3+1 gravity as well. In the case of 
2+1 gravity we demonstrate, that the very existence of black holes makes such 
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2+1 Universe arithmetic. Since most of the Einstein spaces happen to be sym- 
metric, we expect that they are arithmetic in addition in view of the Margulis 
theorem [30]. This possibility is realized in nature only if the black holes exist in 
our Universe. The black holes make our Universe arithmetic. 



2 From train tracks to geodesic laminations 

2.1 Dynamics of train tracks on punctured torus 

Dynamics of pseudo-Anosov homeomorphisms on the four punctured sphere was 
studied in some detail in Ref. [33]. Closely related but more physically interesting 
case is associated with study of self homeomorphisms of the punctured torus. The 
Poincare-Hopf index theorem requires existence of two Y-type singularities, each 
having index -i, as it is explained in Ref. [16]. These singularities can move on 
the surface of the punctured torus thus giving rise to the train tracks dynamics 
as depicted in Fig.l. 

From this picture it follows, that the nontrivial dynamics is effectively caused 
by sequence of meridional r m and longitudinal t; Dehn twists. Using the rules 
set up for dynamics of train tracks [16], one can easily calculate the transition 
matrix by noticing that topologically the state h) is the same as a) while the 
weights on the branches are different. This allows us to write the following 
system of equations 

a = b + 2a, 

b> = c, (2.1) 
c' = a + 2c. 

These results can be neatly presented in the matrix form : 



(2.2) 



The incidence matrix just obtained can be found in Penner's article, Ref. [34], 
where it was presented without derivation. The largest eigenvalue A is found to 
be 

A=i(3 + V5). (2.3) 

Since A > l,this indicates that the dynamics depicted in Fig.l is of pseudo- 
Anosov type. We shall reobtain this result for A below, e.g. see Eq.(2.16), and 
in section 3 using totally different methods. In fact, the presentation above is 
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Figure 1: Fragment of the train track dynamics on once punctured torus 
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only given for the sake of comparison with these new methods to be discussed 
in section 3. These new methods provide the most natural connections between 
the dynamics of 2+1 gravity and the theory of knots and 3-manifolds. 

The reminder of this section is devoted to exposition of some mathematical 
results which will be used in the rest of this paper. 

2.2 The Markov triples 

To begin, let us recall [16,17,19], that a geodesic lamination on a hyperbolic 
surface S is a closed subset of S made of union of disjoint simple geodesies. 
When lifted to the universal cover, i.e. to the Poincare' disc T> whose boundary 
is a circle at infinity, the endpoints of the geodesic lamination determine a 
closed subset (actually, a Mobius strip) 

£ = (S 1 00 xS 1 x -A)/Z 2 (2.4) 

where A is diagonal (x,x),x G 5^ and the factor Z 2 reflects the fact that the 
circle segments representing these geodesies are unoriented, that is the picture 
remains unchanged if the ends of each geodesic which lie on are interchanged. 
This fact has been discussed and used already in our earlier work, Ref. [16] . Since, 
by definition, the lamination is made of disjoint set of geodesies, when lifted to 
T>, there are no circle segments (representing geodesies) which intersect with 
each other. It is intuitively clear, that the dynamics of train tracks should affect 
the dynamics of geodesies. We would like to make this intuitive statement more 
precise. To this purpose, we notice that, when lifted to the universal cover, this 
dynamics causes some homeo(diffeo)morphisms of the circle S^o thus making 
clear the connections with circle maps and mode locking [14,15]. 

If G — tti(S) is the fundamental group of surface S, the endpoint subset £ 
remains invariant under the action of 7r 1 (S l ). Suppose, that the surface S has 
boundaries, e.g. a hole in the case of a torus. Let P C G be the set of peripheral 
elements,e.g. those elements of G which correspond to loops freely homotopic 
to the boundary components. In the case of a punctured torus, G is just a free 
group of two generators: G =< a,b >, and the subset P is determined by the 
commutator aba~ 1 b~ 1 = [a, b] whose trace equals to -2 [35]. Such restriction on 
the trace indicates that the group element [a, b] is parabolic. Parabolic elements 
are always associated with cusps (punctures) on the Riemann surface (for more 
details, e.g. consult Refs.[25,29]).This restriction affects the presentation of the 
group G. Using generators a and b one can construct a "word" W r 

W r = a ai b^ •••a a "6 /3 " (2.5) 

where a\ and (3 r can be any integers while cti and (3 i , i ^ l,r, can be any 
integers, except zero. Not all words thus constructed are different. Those words 
which are conjugate, i.e. W r Wi = WiW r , are considered to be equivalent. If we 
define the equivalence relation as ~ , then the quotient Gj ~ can be identified 
with the set SI of homotopy classes of closed curves on the torus T [35] . Let Cl be 
a subset of SI which corresponds to the nontrivial nonperipheral simple closed 
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curves on T . To account for the presence of a puncture, we require that both a 
and b belong to tt while [a, b] should correspond to the peripheral simple closed 
curve(s).The question now arises: how to find an explicit form of the generators 
a and 6? If we require these generators to act by isometries in the Poincare disc 
T> (or the upper half plane ) we need to find a mapping of the group G into the 
group PSL(2,R). This can be achieved with help of the following identities (in 
the case of SL(2,R)): 

2 + tr[a, b] = (tra) 2 + (trb) 2 + {trabf - tra trb trab (2.6a) 

and 

tra trb = trab + trab^ 1 (2.6b) 

which had been discovered by Fricke [36]. The above identities can be analytically 
extended to SL(2,C) [35]. The projectivized versions of these groups, that is 
PSL(2,R) and PSL(2,C), are groups of isometries of hyperbolic H 2 and H 3 spaces 
respectively. Such an extension is useful for dealing with problems discussed in 
our earlier work, Ref.[25], and will be also discussed later in this work in section 
5. 

Since, as we know already, tr[a,b] — —2, the above identities can be conve- 
niently rewritten as 

x 2 + y 2 + z 2 = xyz (2.7) 

and 

xy = z + w, (2.8) 

where, x — tra, y = trb, z = trab and w — trab . As it is argued by Bowdich in 
Ref.[37] (who, in turn, attributes it to Jorgensen [38]) the first identity, Eq.(2.7), 
is sufficient for restoration of the explicit form of matrices a and b. These are 
given by 

a= 1 -( XZ - y X ) andfe=if yZ X ~ V ). (2.9a) 
z \ x y ) z \ ~y x J 

The above matrices differ slightly from those given in Ref . [37] . This difference is 
essential, however, and originates from the fact that we require deta = detb = 1. 
Such requirement leads automatically to Eq.(2.7) as required. The above choice 
of matrices is not unique since, for example, Eq.(2.7) is also to be satisfied by 
the choice of 

a= 1 -( y X ) (2.9b) 

z \ x xz-y ) 

,ctc. For the integer values of x, y and z the identity, Eq.(2.7), is known as 
equation for the Markov triples discovered in the number theory by Markov [36] 
. Given this equation, one is interested to obtain all integer solutions (triples x, 
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y and z). In mathematics literature, sometimes, related equation is known as an 
equation for the Markov triples [39]. Specifically, one introduces the following 
redefinitions: 

x = 3mi, y = 3m 2 and z — 3m 3 

so that Eq.(2.7) acquires the standard form: 

m 2 + m 2 + m 3 = 3mim 2 m 3 . (2.10) 

The simplest solution of this equation is just mi = m 2 = m 3 = l.To generate 
other solutions it is sufficient to have a "seed" (mi, m 2 , m^) which, by definition, 
obeys the Markov Eq.(2.10). Then, the first generation of Markov triples is given 

by 



(ra'i, 7712,7713), (mi,m' 2 ,m3) and (mi,m 2 ,m 3 ) 



where 



m'i = 3m 2 m 3 — mi; m 2 = 3mim 3 — m 2 ; m 3 = 3mim 2 — m 3 . (2-11) 

Other solutions can be generated via mapping: (mi,m 2 ,m 3 ) — » (— mi,— m 2 ,m3) 
along with cyclical permutation of mi,m 2 and 013. If one starts with (1,1,1) 
then, the Markov tree (its part, of course) is depicted in Fig. 2. 

Emergence of numbers m[, etc., can be easily understood if, instead of 
Eq.(2.10), we would consider the following quadratic form: 

f(x) = x 2 - 3m 2 m 3 x + ml + m 3 (2.12) 

Then, for f(x) = we can identify x = m\ and x' — mj = 3m 2 m 3 — mi. 
In physics literature the dynamics of Markov triples has actually been studied 
already in connection with quasicrystals [24], Id tight binding Schrodinger equa- 
tions (with quasiperiodic potential), etc. [40], and is known as dynamics of trace 
maps. This dynamics can be easily understood based on Eq.s(2.10)-(2.12). In 
short, one studies maps T of the type 

x \ / 3yz — 5 

r- I y -m y I ( 2 - 13 ) 



etc. which possess an integral of motion 

I(x, y, z) = x 2 + y 2 + z 2 - ixyz (2.14) 

invariant under action of T. 

From the theory of Tcichmiillcr spaces [41] it is known, that the length £(7) 
of closed geodesies associated with 7 e G is given by 

tr 2 ( 7 ) = 4cosh 2 (^. (2.15) 
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Figure 2: Fragment of Markov tree close to the original "seed" 
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Let, trj = x (or y, or z), then the Markov triple (1,1,1) corresponds to the 
geodesic whose hyperbolic length is given by 

Z( 7 ) = 2coslT 1 (^) = 21n(i(3 + V5)) = 21nA (2.16) 

where in the last equality use had been made of Eq(2.3). Obtained result pro- 
vides us with the first piece of evidence that the train tracks are directly associ- 
ated with closed geodesic laminations. Since in the case of the punctured torus 
the Teichmiiller space coincides with the Poincare' upper half plane model of 
hyperbolic space, the Teichmiiller distance dr is the same as hyperbolic distance 
Ih [41]. Therefore, we obtain, 

d T = ^\nK = l H = 21nA (2.17) 

which leads to equation 

VK = \ 2 . (2.18) 

That is the stretching factor A is associated directly with the Teichmiiller dilata- 
tion factor K (for more details on K, please, consult our earlier work, Ref. [17]). 

In the light of previous discussion we notice that the spectrum of stretching 
factors is discrete. This is reminiscent already to the energy spectrum of some 
quantum mechanical system. The question arises at this point: is the stretching 
factor A, defined by Eq.(2.3), represents the maximum or the minimum among 
possible stretching factors? A simple calculation based on Eq.(2.15) and Fig. 2 
indicates that A, given by Eq.(2.15), corresponds to the minimum. In the case 
of higher genus Riemann surfaces (perhaps with punctures) Penner [42] had 
demonstrated that the stretching factors are bounded from above and from 
below. 

Let us now obtain the explicit form of matrices a and b, given by Eq.(2.9), 
for the case when x — y = z = 3. An easy calculation produces: 

(2.19) 

yielding the trace of commutator [a, 6] being equal to -2 as required. Introduce 
now two basic matrices 

J J ) and R = ( J 1 ) ( 2 ' 2 °) 

. The letters stand for the "left" (L) and the "right" (R) Dehn twist matrices 
(e.g. see section 3 for more details). Evidently, 

a = RL, (2.21) 

as can be seen by direct calculation. It is also not difficult to check that 

b = L^aR' 1 =L~ 1 RLR- 1 . (2.22) 





/ 2 








-1 




I 







12 



Let us now take a note of the fact that 

L- 1 = ( _\ J ) and JT 1 = ( J - 1 ) . (2.23) 
If we introduce the matrix 

'-±(!~o') <-» 

such that 

*-'-(;!) (2 - 25) 

,then, Eqs.(2.23), (2.24) can be conveniently rewritten as 

L" 1 = IRI (2.26) 

and 

iT 1 = ILI (2.27) 

which also can be checked by direct calculation. Thus obtained relations can be 
equivalently presented as follows: 

LIL = R, RIR = L, RIL = J and LIR=I. (2.28) 

Using the results just obtained, it is clear, that every word W , e.g. see Eq.(2,5), 
can be written in terms of positive powers of L, R and I. Moreover, since 
originally we had only matrices a and b, things can be simplified further. In 
particular, let us consider the word 

W\ = L ai ■ ■ ■ RL ■ ■ ■ L ar R fj r (2.29) 

and the related word 

W 2 = L ai ■ ■ ■ RIL ■ ■ ■ L°' r R l3r . (2.30) 

We can eliminate the contribution RIL in Eq.(2.30) by replacing it with / using 
Eq.s(2.28). If we continue to use Eq.s(2.28) we can evidently get rid of all I 
factors in the middle of the word Wi so that in the end we have to consider only 
the totality of words of the type Wiwith nonnegative integer coefficients. To the 
totality of these words one has to add words like IW\, W\I and IW\I . Since 
in section 4 we shall be interested in the traces of these words, evidently, only 
words of the type W\ and IW\ need to be considered. Incidentally, use of 
Eq.s(2.28) allows us to reduce b, Eq.(2.22), to b = ILRI and, hence, from now 
on we shall use a = RL and b = LR . Such choice is in accord with that known 
in the literature [43,44]. Moreover, to make connections with the results from 
knot theory and 3-manifolds (to be discussed in sections 3 and 5) only words 
of the type W\ should be considered. This peculiarity will be explained further 
below. 
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2.3 The Farey numbers and the Farey tesselation of H 2 



Words of the type W\ are represented by the set of matrices M, 

M=( a J, <«S-7/? = ±1, (2.31) 



7 S 

with integer coefficients. Such matrices belong to the group GL(2, Z) and play 
an important role in the number theory [45]. In the number theory they are 
associated with the Farey numbers. Recall, that the Farey series of numbers $ n is 
the ascending series of irreducible fractions between and 1 whose denominators 
do not exceed n. If p/q and r/s arc two consecutive terms (neighbors) of , 
then 

ps-qr = ±l. (2.32) 

In addition, if p/q , r/s and h/k are three consecutive terms, then the mediant 
r/s is obtained via 

r p± h 

~ = -TT.- 2 - 33 
s q±k 

Hence, — and | in the matrix M correspond to the neighbors in $ n . Given this 
information, one can make two additional steps. First, one can let n go to ±oo 
[46]. Thus, 

#1 : -oo, -1, 0, 1, oo , 



£ 2 : -oo, -2,-1,--, 0, -,l,2,oo , etc. 

Second, by noticing that projectivisation PSL(2,R) of SL(2,R) is just an isom- 
etry of H 2 it is advantageous to map conformally (-00,00) into S^. Once this 
is done, locations of the Farey numbers on S}^ acquire new geometrical mean- 
ing (see also section 5.1). To understand its significance, let us consider the 
Mobius-like transformations associated with matrices L and R, i.e. 

L : z' = ——- and R: z' = z + 1. (2.34) 

z ~\~ _1_ 

Let us notice first that for L-type transformation 

z = 0^z' = 0, z = ±00 — > z' = 1 and z = — 1 — > z' = —00 (2.35a) 
At the same time, for R-type transformation we have 

z = -» z' = 1, z = ±00 -> z' = ±00. (2.35b) 
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Figure 3: Leaky torus in the Poicare disc V model. 



Using these simple results we obtain at once 

6(00) = LR(oo) = 1, o(0) = RL(0) = 1, o(-l) = RL(-l) = -00, 6(-l) = LR(-l) = 

(2.36) 

It is instructive to depict these results graphically, e.g. see Fig. 3. 

Points -00, — 1, 0, 1, +00 on the circle are just the members of 3i .These 
are joined by the circular arcs (not to be mistaken for the hyperbolic geodesies 
which are going to be discussed below). The arrows are in accord with the results 
given by Eq.s(2.36). Identifying sides of the polygon using the arrows depicted 
in Fig. 3 we obtain an orbifold known in physics literature as "leaky torus" [13]. 
To construct the next level of Farey numbers, #2, and, for this matter, etc., 
let us consider the most general modular transformation 

, = az + b ^ 

cz + d w y ' 

where the integer coefficients a, b, c, d are subject to the constraint ad— be = l(to 
be compared with Eq.s(2.31) and (2. 32)). Using this transformation we obtain, 

z = ±00 -> z' = -, z = -» z 1 = -, and z = ±1 -» z' = ±a + b (2.38) 
c d ±c + d v ' 

But the numbers | and | are just the Farey neigbors! Using Eq.(2.33) we obtain 
the following sequence of numbers: 

„ -1-1 1 -1 + 01 + 1 ,„ 1 + 1 

-2= , — = , - = and 2= . 2.39 

+ 1 ' 2 1 + 1 ' 2 1 + 1 + 1 V ; 
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Figure 4: Farey tesselation of V. 



If we place these numbers on the circle we obtain all the members of #2 , etc. as 
it is depicted in Fig. 4. The question arises now: what combinations of L and R 
will lead us to these numbers? Let us consider the generic case of \. With help 
of Eq.s (2.38) and (2.39) we obtain two transformations 

: an d z > = 2i±I , 2a - c = 1 (2.40) 



z + 1 cz + 2 

connecting respectively the numbers 1 and with the Farey number i. Since 
both a and b are integers, the simplest choice is to take a — 1 — c.This then 
produces immediately: 

z' = L(z) and z' = LR(z). (2.41) 

Other examples can now be constructed without problems. Evidently, each 
successive level of can be obtained by some application of combinations of 
L's and i?'s. To make this procedure systematic, let us consider some Farey 
number, say p/q, which has a continued fraction expansion of the type [45] 

-=a H j = [a ,ai, ...,a n ] (2.42) 

q ai + 



1 

a 2 + 



a 3 H 

where a\, a n are some integers. Suppose, we would like to connect the point 
(or 00) with the point p/q. To this purpose, using results given by Eq.(2.23), 
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we need to take into account that 

x 

-x + 1 

Using these results, it is clear, that 

1 ' 



x 1 

x' = LT x x = = and x' = R^x = x - 1. (2.43) 

x 



R- a <>(?-) = L^, (2.44) 

q a i + 7T 

where 

y = 1 —, etc. (2.45) 

1 a 2 H 

Cl 

Using Eq.s(2.43) and (2.44) we obtain, 

L~ 1 R- a o(?.) = j (2.46) 

q oi + — - 1 
h 



and, since 



L~ a x = -j^— (2.47) 

a 

x 



we may as well write, instead of Eq.(2.46), 



iT ai ir ao (-) = — . (2.48) 
q bi 

Taking into account Eq.(2.45) we also get 

L- a2 (L-^Rr^i-)) = a 3 + (2.49) 
V 1 J «i 

Now one has to apply to this i?~ Q3 in order to obtain expression similar to 
Eq.(2.44) so that this process may continue. The end result will depend upon 
wether a n > 1 or a n = 1. When a„ > 1 we can write a„ = (a n — 1) + 1 and then 

use the operator L~( a '> _1 ) to — . Based on the results just obtained, it should 

a n 

be clear by now that it is possible to construct any Farey number starting from 
0,1,-1 or ±oo with the result being some word W of the type 

W = L ai R f3 i ■ ■ ■ L a rRPr (2.50) 

where the exponents ai,f3 1 , ...,a r , f3 r are related directly to numbers m, ...,a n . 
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Consistent triangulation of H 2 (or T> ) invariant with respect to the action of 
PSL(2,Z) can be achieved now. To this purpose first, we connect the point ±00 
with the point j of Fig. 4 with help of the equation 

LRL- 1 R- 1 (-oo) = (2.51) 

which follows easily from Eq.(2.36). Since the points \ and are mediants 
already, it should be clear how to produce the rest of the triangles from these 
two which are basic as it is depicted in Fig. 4. By direct observation, it is easy 
to notice that this figure has an axial symmetry It can be shown [44], that 
the r.h.s corresponds to words of type Wi while the l.h.s. corresponds to words 
of the type IW\ and that these two sets are nonoverlapping. Hence, we shall 
use here and in sections 3 and 4 only the r.h.s. of this figure. Although lines on 
this figure look like geodesies, actually, they are not geodesies in a usual sense 
(so that the statement made in Ref.[47], e.g. see page 566 , that these lines 
are geodesies should be treated with caution ) as it will be explained below. 
Hence, the triangulation in Fig. 4 should be understood only in a topological 
sense (to keep track of how the Farey numbers are related to each other). An 
alternative and very effective geometrical description of Farey numbers (via the 
circle packing in H 2 ) had been proposed by Rademacher [48 ]. We shall discuss 
it briefly in section 5. 

Although the Farey tesselation of T> depicted in Fig. 4 is helpful, it cannot 
be used directly for our calculations since we arc interested in geodesies related 
to the Markov triples. Therefore now we would like to relate these Markov 
geodesies to the Farey tesselation. The simplest Markov geodesies are associated 
with matrices given by Eq.(2.19) (in view if Eq.(2.15)). Since we require these 
matrices to correspond to the closed geodesies on the Riemann surface of leaky 
torus, when lifted to T> (or H 2 ), their ends must lie on (or on the line y=0 
in H 2 model ).The location of the ends of geodesies is determined by the fixed 
points equation which in the case of matrix a is given by 

x = 2 Ji±l = l + ^- T . (2.52) 

X 

Iteration of this equation leads to the periodic continued fraction expansion 
with period 1 [45] 

x = [i] = ^+±. (2.53) 

We had used the notations of Ref . [45] , e.g. see. page 46, to reflect the periodicity. 
Surely, this result can be obtained as well directly from the quadratic equation 

x 2 - x - 1 = (2.54) 

which is equivalent to Eq.(2.52). Eq.(2.54) has two roots 

ari l2 = ^(l±>/5). (2.55) 
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To recover the second root from Eq.(2.52) we have to introduce the following 
change of variable: x=--, in Eq. (2. 52). This then produces after a few trivial 

y 

manipulations 

y = L— . (2.56) 

1 + 



i-y 



Iterating this expression we obtain, 



V = J • (2.57) 

1 + 



1 + - 



This is equivalent to the continued fraction expansion of |(1 — \/5) [45] as 
required. Obtained results are special cases of general theorem proven by Series 
[46] which states that 

x\ = .Too = [711,712, •••] and x 2 = £-00 = (2.58) 

[n ,n_i,n_2, ...J 

so that in both cases we have an infinite continued fractions corresponding 
to irrational numbers. In the case of fixed points for the Markov matrices the 
continued fractions are always periodic [49]. This fact has profound topological 
significance as it will be explained in section 3. In the meantime, we still need to 
clarify the connections between the Farey tesselation of H 2 (or T>) and the results 
just obtained. The comprehensive treatment of this problem by mathematicians, 
surprisingly, had been performed only quite recently [43,47]. Moreover, to our 
knowledge, there are no similar comprehensive treatments for surfaces of higher 
genus (perhaps, with exception of Ref.[50]). The major reason for utilization 
of the Farey triangulation of H 2 is exactly the same as used in approximations 
of irrational numbers by the rationals which belong to the Farey series [45]. 
According to the theory of numbers [45] all rational numbers are equivalent 
since they are connected with each other via modular transformation, Eq.(2.37). 
In our case this means that the vertices of the basic quadrangle depicted in Fig. 3 
upon identification (needed for formation of the leaky torus) all correspond to 
a puncture. Subsequent higher levels of the Farey tesselations do not change 
this result: all Farey numbers correspond to a puncture. Hence, from here it 
also follows that the ends of geodesies which correspond to the nonperipheral 
elements of G are represented by the irrational numbers connected by bi infinite 
sequence of Mobius transformations. In particular case of Eq. (2.52) this becomes 
obvious if we rewrite it in the form 

x = RL{x). (2.59) 

Iterating, we obtain, 

x = RLRLRL ■ ■ ■ (x). (2.60) 
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To clarify the meaning of this result, please, recall that , according to Eq.(2.21), 
RL = a. The eigenvalues X12 of a can be easily found from the equation 

A 2 -3A + 1 = (2.61) 

thus producing 

A lj2 = i(3±V5). (2.62) 

Not surprisingly, Aicoincides with earlier obtained result, Eq.(2.3). These results 
will acquire completely new topological meaning in section 3. 
Consider now the Mobius transformation U\(z) given by 

z' = U x (z)= (^jz = Xz. (2.63) 

This transformation has two fixed points : z* = and z* = 00. Consider yet 
another Mobius transformation W(z) such that 

W(xi) = and W(x 2 ) = 00, e.g. W(z) = z ~ x \ (2.64) 

Z X2 

where x\ and x 2 are given by Eq.(2.55). Hence, using such transformation, we 
obtain, 

WaW-^z) = U x {z), (2.65) 

so that Eq.(2.60) can be rewritten as 

x = WaW^WaW- 1 ■ ■ ■ (x) = U X U X ■ ■ ■ U x (x). (2.66) 

Therefore, indeed, in view of Eq.(2.63), we obtain the expected sequence of 
iterates connecting two fixed points. These points are the initial and the final 
limiting points representing "motion" along the geodesies in H 2 which are just 
semicircles passing through points x\ and X2- It is very important to realize 
that in the case of a torus (and also punctured torus) the hyperbolic H 2 plane 
coincides with the Teichmullcr space. Hence, the "motion" along the geodesic 
in H 2 coincides with the real motion in the Teichmiiller space as was discussed 
qualitatively in our earlier work on 2+1 gravity. More details will be provided 
in sections 3 and 4. 

Consider now the possibility of joining of two Farey numbers by the geodesies. 
That is, let us consider the fixed point of equation 

x= ax + b ad _ cb=1 (26?) 
cx + d 

We need to look for solution of this equation only for the integer values of a, 
b, c and d. It is instructive to discuss some special cases first. Let us begin with 
the case c = b = 0. In this case we obtain, 

x = ^x- ( 2 - 68 ) 



20 



Although this equation looks exactly like Eq.(2.63), these equations are not the 
same since according to Eq.(2.67) we should require ad — 1 which leaves us 
with the only choice: a = d = =Fl.The obtained identity result indicates that 
the hyperbolic-like transformations are not possible if we use just PSL(2,Z). 
Let us now consider the parabolic transformations, e.g. let c = then, given 
that ad — 1, we obtain, 

x = x + b. (2.69) 

This equation has only one fixed point, x* = oo , characteristic for all parabolic 
transformations. Let us now put b = in Eq. (2. 67). Then, again, ad = 1 and we 
obtain, 

(2.70) 



cx + 1 



thus producing another acceptable solution: x* = 0. Finally, let d = or a = 0. 
In the first case we obtain cb=-l so that 



thus leading to the equation 



£1,2 



2 2 

while in the second case we get be = —1 so that 

-1 



ax — 1 

x = (2.71) 



a ± lVa 2 -4, (2.72) 



(2.73) 



x + d 
thus producing 

Xl , 2 = ~±\^d?-A. (2.74) 

In both cases the problem is reduced to finding the Pythagorean numbers, i.e. 
to finding all integer solutions of equation 

a 2 + b 2 = c 2 . (2.75) 

As known results indicate [51], the only solution for d in Eq.(2.74) is d=2 ( 
accordingly, a = 2 in Eq.(2.72)) thus leaving us with just one fixed point. 
Clearly, we can identify thus obtained fixed points with 0, 1,-1 and ±oo as de- 
picted in Fig. 3 since all other points can be obtained by successive applications 
of modular transformations. Hence, the semicircles depicted in Fig. 4 are not 
true geodesies contrary to the statements made in Ref. [47]. Consideration of 
the general case (section 3) produces the same negative result: only one triv- 
ial solution of Eq.(2.75), thus leading to just one (parabolic) fixed point for 
PSL(2,Z). 
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Although wc had provided enough evidence which connects the geodesic lam- 
inations and the train tracks, more systematic treatment of this subject would 
lead us somewhat away from the topics which we had discussed so far. Fortu- 
nately, sections 1.5-1.7 of Ch-r 1 of the book by Penner and Harer [53] contain all 
the required proofs (please, read especially pages 87-101). Much shorter proofs 
could be also found in the unpublished PhD thesis by Lok [54]. Since these re- 
sults are needed only for rigorous mathematical proofs of connections between 
laminations and train tracks, we hope, that our readers will consult these ref- 
erences for better understanding of the obtained results and those which follow 
in the rest of this paper. 

3 From geodesic laminations to 3-manifolds which 
fiber over the circle 

3.1 The Alexander polynomial and surface homeomor- 

phisms 

The Alexander polynomial As(i) for the figure eight knot is known to be [54] 

A 8 (t) = t 2 - 3t + 1 (3.1) 

Consider, quite formally (for the time being only!), the zeroes of this polynomial. 
These are obtained as roots of the equation 

t 2 -3t+l = 0. (3.2) 

A simple calculation produces: 

ti, 2 = ^(3±>/5). (3.3) 

Eq.s (3.2) and (3.3) we would like now to compare with Eq.s(2.61) and (2.62) 
respectively. Since this comparison yields complete coincidence of results the 
rest of this section is devoted to the proof that the above coincidence is not 
accidental. As a result of such proof the connection between the dynamics of 
2+1 gravity and 3-manifolds is naturally established. Unlike Witten's treatment 
of 2+1 gravity [55] which establish this connection through reformulation of 
this problem in terms of the Chern-Simons field theory, our treatment does not 
require field-theoretic arguments at all and is based mainly on Thurston's theory 
of 3-manifolds [19]. The condensed summary of relevant results of Thurston [19] 
and McMullen [56] is given in our earlier publications, Ref.[16,17]. 

We would like to begin with reviewing some facts about the Alexander poly- 
nomial. Although our earlier published review [57] provides sufficient physical 
background on knot polynomials and, in particular, on the Alexander polyno- 
mial, this background is not sufficient for our current purposes. Hence, we shall 
avoid references to physics literature when we shall talk about the Alexander 
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polynomial. To simplify matters, we shall treat only the case of the Alexan- 
der polynomials for knots. The case of links is considerably more complicated 
and will be treated in a separate publication. For the case of knots, let V be 
the Seifert matrix of linking coefficients [18,58], then the Alexander polynomial 
Ax(t) for a knot K is given by 

Ax(t) — dct(V T — tV) (3.4) 

where V T is the matrix transpose of V. Such defined polynomial has some 
additional remarkable properties [18,58]. For instance, 

A K (1) = ±1 (3.5) 



AkW+AkH- 1 ) (3.6) 

where the symbol = denotes an equality up to a constant multiplier (e.g. see 
Eq.(3.1) for an obvious example). The above properties of Ak allow us to write 
it as polynomial of even degree 2r with integer coefficients dj : 

1r 

Ax{i) = 51 cii? , a 2r -i = a t . (3.7) 
»=o 

Using the inversion symmetry property reflected in Eq.(3.6) the following Lau- 
rent expansion for Ak can be written: 

A K (t)=a r + a r+1 (t + tr 1 ) + ... + a 2r (t r + tr r ). (3.8) 

We are not interested in all possible knots and their Alexander polynomials since 
not all knots are of relevance to surface dynamics. As we had discussed earlier 
[17], only fibered knots are of relevance. The easiest way to talk about fibered 
knots embedded in S 3 is trough consideration of their complements S 3 \ K in 
S 3 . These are just 3 manifolds fibering over the circle S 1 . Indeed, let S be 
the Seifert surface associated with knot K. Incidentally, the punctured torus 
is the Seifert surface for both the figure eight and the trefoil (right and left) 
knots [18,58]. There is another knot, the bridge knot b(7,3), which also has the 
same Seifert surface but it cannot be fibered over the circle [59] . The Alexander 
polynomial At for the trefoil is known to be [54] 

A T (t) = t 2 -t + 1. (3.9) 

The zeros of this polynomial are readily obtained: 

ti,2 = ^(l±n/3). (3.10) 

Apparently, they have nothing to do with the discussion made in the previous 
section. This, however, is not true as we shall soon demonstrate. To this purpose 
consider an orientation preserving surface homeomorphism h: S — > S. In the 
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case of the punctured(holed) torus the homeomorphism should respect the pres- 
ence of a hole. The circumference of this hole is just our base space S 1 (which, 
as it is not too difficult to guess, is just our knot K since the knot is just a 
circle embedded into S 3 ). The Seifert surface itself is a fiber and the 3-manifold 
is just a fiber bundle constructed in a following way. Begin with products 5x0 
(the initial state) and Si x 1 (the final state) so that for each point x € S we 
have (x, 0) and (h(x), 1) respectively. The interval / = (0, 1) can be now closed 
(to form a circle S 1 ) by identifying with 1 which causes identification: 



The fiber bundle (also known in the literature as mapping torus [56,60]) 



is the 3-manifold which fibers over the circle and is complementary to the fibered 
knot in S 3 . The interval (0, 1) can be associated with some local time. The cyclic 
character of the process leading to formation of 3-manifold(s) is not essential 
as it will be demonstrated below. Therefore, actually, the time interval can be 
taken from -co to oo . The periodicity naturally occurs, if homeomorphisms are 
associated with motion along the Markov geodesies in the Teichmiiller space, 
e.g. see Eq.(2.60). Indeed, the homeomorphisms of surfaces are associated 
with dynamics of train tracks. If we start with train track dynamics, it is 
in one-to one correspondence with dynamics of geodesic laminations and this 
dynamics, in turn, is associated with motions in Teichmiiller space, e.g. along 
some geodesies in this space. Hence, the periodicity occurs quite naturally. 
Othal [61] and McMiillen [56] had proved that the situation just described for the 
punctured torus persist for Riemann surface (with marking and /or boundaries) 
of any genus g as it was briefly mentioned earlier in our work, Ref.[17]. More 
specifically, they proved the following theorem 

Theorem 3.1. Let ip : S — > S be a pseudo-Anosov homeomorphism of com- 
pact surface with negative Euler characteristic. Then, in order for the mapping 
torus T^f, to have a hyperbolic structure one is looking for the related hyperbolic 
manifold = S x R on which the homotopy class of ip is represented by an 
isometry a . Then, M^j < a > is homeomorphic to T^. 

It can be demonstrated, e.g. Proposition 5.10 and the comments which follow 
in Rcf.[18], that classification of all fibered knot complements can be formulated 
in terms of fibering surfaces and maps of such surfaces. Therefore, naturally, this 
fact is reflected in the associated Alexander polynomials. Unfortunately, the 
conjecture that all fibered knots are classified one-to-one by their Alexander 
polynomials happens to be wrong. Morton had demonstrated [62] that for the 
Seifert surfaces of genus g>l there are infinitely many different fibered 
knots for each Alexander polynomial of degree >2. This fact by no means 
diminishes the role of the Alexander polynomial in dynamics of 2+1 gravity, 
it just makes knot interpretations [22,55] of gravity less convincing. Similar 



(x,0) = (h(x),l). 



(3.11) 



T h = (SxI) /h 



(3.12) 



24 



negative conclusions had been reached recently with help of absolutely different 
set of arguments in Ref.[63]. 

From the knot theory [18 ] it is known, that instead of expansion, Eq.(3.7), 
for the Alexander polynomial for any knot, in the case of fibered knots one 
has to use 

& K (t) = j2 a * f > ( 3 - 13 ) 

i=0 

where g is the genus of the associated Seifert surface. Moreover, 

Ajf(O) = a = a 2g = ±1, (3.14) 

,i.e. the Alexander polynomial is monic. Using Eq.(3.4) we obtain at once 

A K (0) = det(U T ) = det(U) = ±1 (3.15) 

to be compared with determinant in Eq.(2.31). Using Eqs.(3.15) and (3.4) we 
obtain without delay, 

A K (t) = dctiV^V 71 - tE) = det(M - tE), (3.16) 

where E is the unit matrix and M is the monodromy matrix responsible for 
surface automorphisms [18, 58]. This fact easily follows from the observation 
that A K (0) = det(M) = ±1. 

Remark 3.2. Equating the Alexander polynomial to zero and finding the 
roots of this polynomial is equivalent to solving the eigenvalue problem for 
surface automorphisms which produces the stretching factors allowing one to 
distinguish between the hyperbolic, pseudo-Anosov-type ( when the largest root 
is real and greater than one), and the periodic, Seifert fibered-type (when the 
modulus of the largest root is equal to one), regimes of surface homcomorphisms. 

Let us discuss Eq.(3.16) a bit more. It can be easily shown [64] that, if 
M e GL(2,Z),then 

A K (t) = t 2 - (trM)t + dctM. (3.17) 

But, as we know already, detM = ±1 . Hence, in this case we obtain, instead 
of Eq.(3.17), 

A K (1) =t 2 - (trM)t± 1. (3.18) 

Using Eq(3.5) we obtain as well 

Ajf(l) = 1 - {trM) ± 1 = ±1. (3.19) 

This leaves us with two options: tr(M) = 1 or 3. In the first case we reobtain 
the Alexander polynomial for the trefoil, Eq.(3.9), and in the second for the 
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figure eight, Eq.(3.1), knots. No other options are available! Hence we had 
just proved (in a somewhat different way as compared with Ref.[18 ]), that the 
trefoil and the figure eight knots are the only two fibered knots associated with 
the Sefert surfaces of genus 1. Moreover, the conditions trM = 3 and detM = 1 
lead us directly to the Markov matrix a, Eq.(2.19). This matrix, Eq.(2.60), and 
the Theorem 3.1 provide direct connection between the hyperbolic 3-manifold 
associated with complement of figure eight knot and the pseudo-Anosov surface 
homeomorphisms associated with a. These conclusions are in accord with results 
of Thurston [19], e.g see section 4.37, where they had been obtained in a different 
way. 

At the same time, the conditions trM — 1 and detM = 1 are also very 
interesting since they are associated with the matrix 

M =(-i J) OT (1 ( 3 - 2 °) 

whose eigenvalues we had calculated already in Eq. (3.10) . It is easy to check that 
such surface homeomorphism is not associated with motion along the hyperbolic 
geodesic since the fixed point equation 

x = 1 - - (3.21) 

x 

which is equivalent to Eq.(3.9) (if we require Ax(t) to be zero) does not have real 
roots. Using the same methods as in section 2 we easily obtain that projectively 
the matrix M is equivalent to the combination of LR^ 1 . Such transformations 
do not fit the Theorem 3.1 and the associated 3 manifolds are known as Seifert 
fibered spaces. Since the trefoil knot is the simplest representative of torus 
knots, it can be demonstrated [18], that the complements of all torus knots 
in S 3 are associated with the Seifert fibered spaces. According to Thurston's 
classification of surface homeomorphisms those leading to the Seifert-fibered 
spaces are known as periodic. The detailed structure of the periodic phase 
may be very complicated [65-67] and, hence, very interesting from the point 
of view of physical applications. The discussion of all emerging possibilities, 
surely, requires separate publications. The geometrical richness of the Seifert- 
fibered regime provides a likely explanation of different states of order in the 
case of liquid crystals as we had briefly mentioned in the Introduction and earlier 
works. E.g. solid and hexatic phases might be associated with the Seifert-like 
while liquid and/or gas phases may be associated with the pseudo-Anosov- like. 
In this paper we shall be concerned only with transition between the pseudo- 
Anosov and the Seifert-like phases leaving the detailed analysis of possibilities 
emerging in the Seifert-fibered (periodic) phase for future publications. In order 
for our results to be consistent with the rest of the literature on phase transitions 
we need now to introduce several new concepts. 
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3.2 Mahler measures and topological entropies 

Following Ref.[68] let us consider a monic polynomial with integer coefficients 

F(x) =x d + a^x 4 ' 1 + ... + a lX + a (3.22) 

where ao = ±1. If on are the roots of this polynomial, then, equivalently, we 
can rewrite it as 

d 

F{x) = Ylix-on). (3.23) 

i=l 

The logarithmic Mahler measure m(F) can be defined now as 

m(F)= In M(F) (3.24) 

where M{F) is given by 

d 

M(F) = JJmax{l,|aj|}. (3.25) 

i=l 

The following theorem is attributed to Walters, Ref.[69], section 8.4. 

Theorem 3.3. The topological entropy of the transformation M is equal to 
the logarithmic Mahler measure of the characteristic polynomial of the matrix 
M. 

To make a connection with statistical mechanics, we would like to notice 
that the fact the polynomial F(x) is monic is not at all restrictive (we used it 
only in order to make connection with the previous discussion). Let us recall 
at this point some facts from Lee and Yang theory of phase transitions [7] . For 
instance, for a gas of 9Jt atoms which can be packed into the volume the grand 
canonical partition function 5 is given in a usual manner as 

sot 

£ = £^". (3.26) 

n=0 

Hence, the grand partition function is just some polynomial in fugacity z. Surely, 
the polynomial, Eq.(3.26), must have some zeros. These zeros can be only in 
the complex plane C and they have to come in pairs of complex conjugates for 
Eq.(3.26) to be real for real z's. Hence, effectively, we can rewrite H as follows 

o m/z 

n— 1 

Since 5(2 = 0) = l,the above result can be conveniently rewritten as 

SOT 

Z(z) = H(l--)(l-^). (3.28) 



71=1 
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For physical applications one is interested in finding of Zns that is 

m/z 

$(z)=]nE= Hz 2 - 2cos0„ + 1) (3.28) 

n=l 

where use had been made of the fact that the complex zeros should lie on the 
unit circle in the complex plane C (this is demonstrated below). If g(0) is the 
density of these zeros, one can replace the above summation in Eq.(3.28) by the 
integration with the result 

${z)= j d6g{6) In | (z 2 - 2 cos + 1) | (3.29) 
o 

which is written with account of symmetry of the integrand. To make a connec- 
tion between this result and the logarithmic Mahler measure it is only sufficient 
to use 

Lemma 3.4 (Mahler's lemma) [68] For any nonzero F <G C[x], 

l 

m(F) = J d0m|F(cxp(27rz0))| . (3.30) 
o 

Proof. This result follows at once if one can prove Jensen's formula: for 
any a € C 

l 

J ln|exp(27n0) - a\ = lnmax{l, |a|}. (3.31) 
o 

The detailed proof of this result can be found in Ref. [68]. Using Eq.s(3.23)-(3.25) 
the Lemma 3.4. is proved. 

From here, the following theorem follows at once. 

Theorem 3.5. m(F) = if and only if \ai\ = 1 for all i in Eq.(3.23), 
provided that F{x) is monic polynomial. 



This result coincides with that of Lee and Yang, e.g. see Theorem 3, e.g. 
see page 418, Eq.(57) of Ref. [7], where it was obtained using different methods. 
Theorem 3.5. does not imply that for all monic polynomials m(F) = O.This 
is definitely not the case for the Alexander polynomial, Eq.(3.1), for the figure 
eight knot. Evidently, for such knot the topological entropy is given by 

m 8 = lni(3 + \/5) (3.32) 
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while for the trefoil knot m-r, is, indeed, zero in view of Eq.(3.10). Hence, the 
logarithmic Mahler measure can be used for the Alexander polyno- 
mial of any fibered knot to distinguish between the pseudo-Anosov 
and the periodic regimes. We would like now to complicate matters by con- 
sidering 

3.3 The incompressible surfaces in the once punctured 
torus bundles over S 1 

Although the previous analysis is interesting in its own right, the situations 
which we had considered so far (with the trefoil and the figure eight knots) are 
not the only possibilities for the punctured torus automorphims. To go beyond 
these possibilities means to abandon almost completely the connections with 
knots and to concentrate attention on solutions of the equation 

D (i) = det(M - tE) = 0. (3.33) 

Since we now do not require 2) (t) to be connected with knot polynomials, there 
is no need for D (1)=±1. At the same time, Eq.(3.33) is legitimate equation for 
finding of stretching factors characteristic for pseudo-Anosov homeomorphisms. 
For an arbitrary matrix M e SL(2,Z) we obtain, instead of Eq.(3.18), the 
equation ( n € Z ) 

t 2 - nt ± 1 = (3.34) 
which produces the following roots: 



ti, 2 = -(n±y/n 2 T4)- (3-35) 

From here, we see that only n which obey \n\ < 2 produce the stretching factors 
characteristic for the periodic (Seifert fibered) phase. For all |n| > 3 the stretch- 
ing factors are characteristic of the pseudo-Anosov phase. These results can be 
interpreted topologically with help of the notion of incompressible branched 
surfaces introduced by Oertel [70]. We would like to discuss these surfaces in 
order to make connections with the results of section 2 and in order to provide 
topological interpretation of the results of section 4. 

The notion of branched surfaces is inseparably connected with the notion of 
train tracks. As soon as we would like to visualize the dynamics of train tracks 
in time, we encounter branched surfaces which are perpendicular to the original 
surface as it is depicted in Fig. 5. 

Naturally, these branched surfaces inherit the weights, e.g. see Fig.l, from 
the associated with them train tracks and these weights obey the switch condi- 
tions as for the train tracks. This is depicted in Fig. 6a). The weights are some 
positive integers which could be interpreted as number of surfaces in the stack, 
Fig. 6b). These branched surfaces inherit from the train tracks some topological 
restrictions as depicted in Fig. 7 
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Figure 5: Branched surfaces associated with dynamics of train tracks 



Figure 6: a) Kirkhoff -type rules for branched surfaces . b) Another interpreta- 
tion of the weights on branches (read the text for details) 
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Figure 7: Topological restrictions which branched surfaces inherit from train 
tracks. 

For instance, in the case of train tracks the monogons are forbidden as well 
as discs of contact, etc. (for details, please, contact Ref.[53]). Floyd and Ortel 
[71] (see also Oertel [70 ] , page 387) had proved the following 

Theorem 3. 5. There is a finite collection of incompressible branched sur- 
faces in some 3-manifold M such that every 2-sided incompressible, d -incompressible 
surface in M is carried with positive weights by a branched surface 03 of the col- 
lection. 

The surface S is called 2-sided in M if there is a regular neighborhood of S 
homeomorphic to S x I where J=[0,1], [72]. If S ^ S 2 , RP 2 , or a disc D 2 which 
can be pushed to dM, then £ is incompressible if every loop (a simple closed 
curve) on S which bounds an (open) disc in M \ S also bounds a disc in S. If 
S has boundary, then S is also d — incompressible if every arc a in S (with 
d(a) <G dS) which is homotopic to dM is also homotopic in S to dS [19], section 
4.10. 

Fig. 6b) shows the local model for the fibered regular neighborhood 7V(*B) of 
*B , it also depicts the horizontal boundary <9/jiV(Q3) and the vertical boundary 
d v N(<B). 

The branched surface 25 is incompressible if : 

a) *B has no discs of contact; 

b) 9/^(03) is incompressible and d— incompressible in M\N where TV is 
just an interior of N(^8) and d— compressing disc for <9/j./V(Q3) is assumed 
to have boundary in dM U 9^iV(Q3); 

c) there are no monogons in M\N . 

Being armed with such definitions, we would like to show how actually these 
surfaces can be built. This task was actually accomplished in Refs. [73,74 ]. Al- 
ternative treatment can be found in recently published paper by Minsky [47]. 
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It is essential to contact these references for detailed study of this very involved 
topic. The theory of incompressible surfaces in connection with general theory 
of 3- manifolds can be found in the review by Jaco [75]. Here, we restrict ourself 
only with very basic facts which are needed for topological interpretation of the 
results which follow in section 4. 

The mapping torus construction, Eq.(3.12), does produce 3- manifolds which 
are the complements of the fig. 8 and trefoil knots respectively. It is clear, based 
on Eq.(3.19), that nothing much beyond this can be obtained. To get more, 
one needs to generalize this construction. Following, Ref.[73], we would like to 
consider more involved way of constructing torus bundles. For example, let S 
be some surface and let X\ , X 2 be the topological spaces associated with the 
mapping n / i : S x [0, 1] — > , i=l,2. For some surface homeomorphism h: S — > 
S define Xi jh^-i to be the quotient space obtained from X\ U^2 by identifying 
7 1 (x, 1) with 72 0). In this picture the mapping torus Th, Eq.(3.12), is just 
a quotient space obtained from X\ by identification of r y 1 (x, 1) with -f 1 (h(x), 0). 
Now, however, we can extend this construction by letting i to range from 1 to 
n.Thus, if 7j : S x [0, 1] — ► Xf and hi : S — > S are homeomorphisms for i=l,..., 
n, then the fiber bundle 

X = Xi/ h 1 X 2 / h 2 • • • /h n - 1 X n /h n (3.36) 

is associated with some 3-manifold which, actually, can be obtained from the 
figure 8 3-manifold by means of hyperbolic Dehn surgery [19,76]. Before ex- 
plaining its meaning, several simpler concepts need to be elucidated. First, 
since according to section 2, all toral homeomorphisms can be performed with 
help of the right R and the left L Dehn twists (for orientation preserving home- 
omorphisms), the fiber bundle construction just described is equivalent to some 
transformation ip of the punctured torus given by 



ip = R ai L a2 ■ ■ ■ L a2n , m>l, (3.37) 

in accord with Ref. [73] . Since each torus is completely determined by the ratio 
r of its sides, e.g. r, = a%jbi , it is clear, in view of Eq.(3.36) that the transfor- 
mation ip should be such that ip(a,i/bi) = di+2n/&j+2n for all i and fixed n This 
prompts us to discuss the properties of the partition function Z = trip which we 
postpone till section 4. In this section we would like to discuss a bit more the 
topological issues. According to the results of section 2 the transformation ip is 
of the same kind as the transformation W in Eq.(2.50). This means that the cu- 
mulative result of matrix multiplications in Eq.(3.37) is just some matrix M of 
the type given by Eq.(2.31). In section 2 we had already discussed special cases 
of equations for geodesies. Here, we would like to investigate the most general 
case. Using Eq.(2.67) we obtain the following roots which provide locations of 
the begining and the end of hyperbolic geodesies on 

xi,2 = ^-(a-d±^{a + d) 2 - 4). (3.38) 
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Since a + d is just the trace of M we have to impose the usual requirement that 
tr 2 M > 4 for the transformation to be hyperbolic. In addition, however, we have 
to take into account that tr 2 M is some nonnegative integer. As we had argued 
in section 2, the roots Xi t 2 cannot be rational numbers. Hence, whatever they 
might be, they will belong to some quadratic irrationalities. From the number 
theory [45 ] it is well known that the continued fraction expansions of quadratic 
irrationalities are periodic as it had been mentioned already after Eq.(2.58) and 
this fact also explains the periodicity of ip transformation. 

Remark 3.6. Since, according to section 2, the exponents ai,...,a 2n in 
Eq.(3.37) are in one to one correspondence with the coefficients of continued 
fraction expansion, this periodicity is in one to one correspondence with the 
fiber bundles constructed in Eq.(3.36). 

That this is indeed the case is shown in great detail in Ref. [73,74]. More- 
over, since the transformation (p is hyperbolic, it acts by translations along the 
geodesies in T> whose ends are determined by the roots of Eq.(3.38). This has 
been illustrated already in Eq.s (2.60) and (2.66). Imagine now that we move 
along one of such geodesies, say (5, that is we move in the Teichmiiller space 
of the punctured torus. Such motion will be associated with the sequence 
of crossings of triangles of the Farey tesselation. To understand why this is so 
the folowing arguments are helpful. Each torus can be triangulated that is it 
can be represented as at least two triangles. Obviously, one of this triangles is 
sufficient for complete characterization of the torus. The triangle is determined 
by the slopes of its sides. E.g. the undistorted triangle is determined by the 
triple A(T) =< 1/0, 0/1, 1/1 > . It is sufficient to take a look at Fig.s 3 and 4 in 
order to recognize one to one correspondence between this triple and the Farey 
numbers 3i for the largest triangle in the Farey tesselation. Hence, motion in 
the Teichmiiller space of the punctured torus is indeed associated with sequence 
of triangles in the Farey tesselation of T> since M(A(T)) = A(M(T)) where 
M{z) is defined by Eq.(2.37). Topologically, we can visualize this chain of tri- 
angles as an infinite strip , e.g. see Fig. 8, which consists of triangles crossed 
by ©.The numbers aj > 1 indicate the number of smaller triangles within each 
larger triangle. They are in one to one correspondence with the exponents in 
Eq.(3.37). The periodicity of the transformation ip is reflected in the periodicity 
of the triangulation pattern in the strip J2 V -To make all these statements more 
physical, we need to discuss the meaning of hyperbolic Dehn surgery (actually, 
Dehn filling!) now. 

Let M be a hyperbolic 3-manifold whose boundary DM is a torus T. If we 
write T=R 2 /Z 2 and choose for basis vectors ei,e 2 of the lattice Z 2 then, for 
any coprime pair of integers (p,q) the element pei + qe 2 determines some simple 
curve on T. Consider another solid torus S 1 x D 2 . This torus we glue into M 
in such a way that the curve pei + qe 2 in dM is being glued to the meridian 
(that is to the curve of slope 1/0) of the solid torus. 

Remark 3.7. It is being said that this new 3-manifold M(p,q) is obtained 
from the old M by the operation of Dehn filling along the curve of slope p/q. 
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Figure 8: Bi infinite strip of edge-paths associated with hyperbolic trans- 
formation (p 

Hatcher has proved the following remarkable theorem (using branched sur- 
faces) [77]. 

Theorem 3.8. A 3-manifold with a single torus boundary component has 
only finitely many boundary slopes. 

Remark 3.9 a). From the previous discussion it follows that there should 
be a correspondence between the incompressible surfaces and the boundary 
slopes. b) The boundary slopes are also in correspondence with QU{oo} and, 
hence, with the Farey tesselation of V. 

This correspondence can be rephrased in physically familiar terms as we 
would like to demonstrate now. To this purpose, let us recall that the continued 
fraction expansion for a which is quadratic irrational can be written as 

a = [do, a r , a r+ i, a r+ 2 n , »r+i, ■■■a r +2n, ■■■]■ (3.39) 

The quadratic irrational is purely periodic [78 ] if we have a m +2n = »m for all 
m (i.e. in this case, the terms ao, a. r are absent ). Evidently, Fig. 8 represents 
just this case. Let us now associate the boundary slope p/q with continued 
fraction 

p/q=[a 1 ,...,a 2n ] (3.40) 

(e.g. see the discussion after Eq.(3.37) ) and consider directed random walks 
from the point 1/0 to the point p/q on the diagram, Fig. 8. The random walk 7 
is directed if it does not have backtracks, moreover, the walk is minimal if two 
successive steps never belong to the same triangle. A brief introduction to the 
properties of such walks can be found in our earlier work, Rcf. [79] , where such 
walk was used for description of the discretized version of the Dirac propagator. 
Such discretization is useful in some problems relevant to physics of semiflexible 
polymers. 
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Following our earlier work, Ref. [79] , let N=n be the total number of steps in 
the directed(edge-paths) walk 7 and let N + be the number of right turns while 
N_ be the number of left turns in such a walk. Associate a system of N Ising 
spins U{ with the walk so that Oi = +1(— 1) will correspond to the turn to the 
left (right). With such defined rule for spins, the "magnetization" M is given 

by 

N 

M = £Vi. (3.41) 

i=i 

In the case if eigenvalues of the matrix ip , Eq.(3.37), are positive the boundary 
slope m(n) is given by [74,80] 

m(n) = ^M. (3.42) 

If the matrix <p has negative eigenvalues then, instead of Eq.(3.42), one has 
m(n)=(l/4)M+l/2. In the next section thus introduced quantities acquire new 
statistical mechanical meaning. Floyd and Hatcher [74] have proved the fol- 
lowing theorem 

Theorem 3.10. If ip (Eq.3.37) is hyperbolic, then a connected, orientable 
incompressible, d— incompressible surface in M v is exactly one of: 

a) the peripheral torus dM v , 

b) the fiber T 2 - {0}, 

c) a finite number (>2) of non-closed surfaces S 7 indexed by minimal 
edge- paths 7. 

From here, it follows, in particular, that for transformation given by the 
matrix a, Eq.(2.21), which is relevant to 3-manifold associated with the com- 
plement of fig. 8 knot, there are exactly 2 non-closed orientable incompressible 
d— incompressible surfaces. According to Thurston (Ref. [19], section 4.40) these 
are just T 2 \2disks surfaces which are oppositely oriented. The above Theorem 
3.10 does not include the nonorientable incompressible surfaces. This case was 
studied by Przytycki [80,81] whose work extends the work by Floyd and Hatcher 
and is based on the same edge-path methods. For the complement of fig. 8 knot 
3-manifold the incompressible unoricntcd surface is just the Klein bottle-disk as 
it was also shown by Thurston in the same section 4.40. We refer our readers 
to the original papers for more details. In the meantime, we would like now to 
provide some physical interpretation of the obtained results using some recent 
results from statistical mechanics of number-theoretic spin chains. 

4 Thermodynamics of the Farey spin chains and 
statistical mechanics of 2+1 gravity 

In our previous work, Ref. [16], we have used meanders and Peierls-type argu- 
ments for description of phase transitions in liquid crystals and gravity. Such 
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type of approach, although provides some estimate of transition parameters, is 
not well suited for more refined analysis since even the notion of the order pa- 
rameter, central to all theories of phase transitions, is not so easy to implement 
within such an approach. In the previous section we had introduced the log- 
arithmic Mahler measure, Eq. (3.24), which is ideally suitable for description 
of dynamical phase transitions and in, addition, perfectly fits Yang and Lee 
theory of phase transitions [7]. To make use of this measure in this section, we 
need to provide an equivalent definition of this measure now (incidentally, this 
alternative definition provides an extra link between the statistical mechanics 
and the number theory [68]). To this purpose, instead of Eq.(3.23), introduce 

d 

A„(F) = n«-l)- (4-1) 

1=1 

Using thus defined polynomial A„ it can be shown that, provided no zero of 
F(x), Eq.(3.22), is a root of unity, 

m(F) = lim -ln|A„( J F)| . (4.2) 

n^oo TL 

In the case if polynomial F(x) has zeros which are roots of unity, m(F) — 
0. Being armed with such results, following Refs[82,83], let us introduce the 
spin-like variable <r, = {0, 1} , i=l-k , which is related to the the Ising spin 
Si=(-l) ai . Define now inductively the matrix Dk via the following set of rules: 

a>=(; ;) (4.3) 

and 

Dk = L 1 ~ ak R ak D k _i(a\, <Tfe_i). (4.4) 
The energy Ek can be defined now as 

S fc -lnr fe (4.5) 
with Tfc = Tr(Dk). This allows us to introduce the partition function 

Zkifi) = Yl exp(-/3£ fe ). (4.6) 
Wi} 

with control parameter (3 playing a role of the inverse temperature. The free 
energy 3^/3) can be defined now in a usual way as 

= lim F k (4.7) 

K—>00 



where 



F fe = -^ln(Z fe (/3)). (4.? 
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To make connection with Eq.(4.2) we notice that 



m(F) = lim \e u . (4.9) 

k— >oo K 

Since E k is random variable, more appropriate quantity is the average Mahler 
logarithmic measure which can be interpreted as an average energy U(f3), i.e. 
we have 

U((3) = lim f4(/?) with U k ((3) = JL(/3F k (/3)). (4.10) 

fc— »oo Cp 

The average magnetization (per site) associated with the average boundary 
slope, e.g. see Eqs.(3.41),(3.42), can be defined accordingly as 

k 

M k(P)=<lj2^>k (4.11) 

i=i 

where, as usual, for any observable O the average < ... >fcis defined by 

<0> k 08) - ^= . ~z - . (4.12) 

J2 cxp(-/3£ fc ) 

W\ 

Detailed calculations performed in Refs.[7-10, 82, 83] allow us to avoid repe- 
titions. Hence, we only provide here the summary of results. In the thermo- 
dynamic limit, co, the partition function, Eq.(4.6) acquires the following 
form: 

announced earlier in the Introduction, e.g. see Eq. (1.4). 

Remark 4.1. The partition function Z(f3) can be brought into form which 
coincides with that for lattice gases, e.g. see Eq.s(3.27) and(3.28). To this 
purpose let us introduce the notation = T(^)(f3— l)ir~ z((f3) where T(x) is 
just the Euler's gamma function. Then, following Riemann [2], we notice that 

p 

where p ranges over the roots of the equation £(/?) = 0. Hence, the complete sta- 
tistical mechanics treatment of 2+1 gravity depends crucially on our knowledge 
of distribution of zeros of Riemann zeta function g(p) and, hence, on construc- 
tive solution of Riemann hypothesis. 
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Remark 4.2. Taking into account that the Riemann's Zeta function has 
the following integral presentation [1,2] 



C(/3) = ivm / dx 



r(/3) J e x -l 
o 

one can interpret the obtained results in terms of the thermodynamic properties 
of some fictitious Bose gas so that the phase transition in such gas to some extent 
resembles Bose condensation. This topic is discussed further in the Appendix 
A.l. 

Since solution of the Riemann hypothesis is not yet found (see, however, 
Ref.[5]) we may be content for now by finding some meaningful estimates. For 
instance, Eq. (4.13) allows one to obtain the critical "temperature" (3 cr = 2 
[ 9 ]. In the high temperature (pseudo-Anosov) phase, i.e. for < (3 < (3 cr 
, the averaged logarithmic Mahler measure U{(3) is bounded by the following 
inequalities [9]: 

In 2-3 In § 3 

2 _p 2 <CT(/?)<ln-. (4.14) 

Since for 1.7 <(3 < 2 the l.h.s. of this inequality is negative, more accurate 
estimate was obtained for 1< (3 < 2 : 

U(0)> \iPcr-0)- (4-15) 

With such improved estimate Contucci and Knauff [9] had demonstrated that 
U(P) > for 0< (3 < 2. Thus, the high temperature phase is indeed of pseudo- 
Anosov type. For (3 > (3 cr , i.e. in the low temperature (periodic or hexatic 
[20] ) frozen phase, U(f3) — 0, which is physically meanigful. In this phase also 
3"(/3) = 0. At the same time, the average magnetization 

M{(3) = lim M k (4.16) 

k — >oo 

is equal to one in the low temperature phase and is zero in the high temper- 
ature phase. This also makes physical sense since, according to the results of 
Przytycki, e.g. see Proposition 3.2. of Ref. [81], new M(p,q) manifolds ob- 
tained by Dehn filling from 3-manifold M which is complement of the fig. 8 knot 
are Seifert fibered only if p/q=±l/l, ±1/2, or ±1/3. In order to connect this 
result with physics it is important to realize that the "magnetization" defined 
in Eq.(3.41) is made of spins which are periodically arranged in the bi infi- 
nite strip ^2 . The thermodynamic averages are appropriate only if the period 
n^ oo. In the case of finite n some n-dependent corrections are usually ex- 
pected [20] for all observables. These corrections naturally should vanish in the 
thermodynamic limit n— > oo. Hence, identification of the average magnetization 
with the average slope cannot be done automatically. This casts some doubts 
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Figure 9: Graph ^ dual to the Farey tesselation of V. 



on the appropriateness of the above thermodynamic formalism for description 
of statistical mechanics of 2+1 gravity. Fortunately there are other ways to 
arrive at the same conclusions which we are now going to discuss in this and 
the following section. 

First, instead of considering the Farey triangulation of T> we can consider 
the associated with it dual tree as depicted in Fig. 9. Following Bowditch [35 ], 
we notice that, actually, topologically such tree is in one to one correspondence 
with the Markov tree depicted in Fig. 2. This observation can be used for de- 
scription of random walks on such trees. This walk can be associated naturally 
with motions in the Teichmuller space of the punctured torus as has been ob- 
served already by Penner [84]. Let ^2 denote the dual graph depicted in Fig. 9 
and let VQ^), EQ^) and denote respectively the sets of vertices, edges and 
the regions complementary to ■ Each vertex lies at the boundary of three 
complementary regions X, Y, Z G while each edge meets four complementary 
regions X, Y, Z, W. The subset f2 of nontrivial nonpcripheral closed curves on 
T introduced in section 2 can now be identified with f2 because X,Y, Z,W G O 
correspond respectively to the equivalence classes of generators a,b,ab and ab^ 1 
of the free group G as discussed in section 2. Looking at the Fig. 9 it is easy to 
recognize that the regions of Q are in one to one correspondence with rationals 
QU{oo} of the Farey tesselation. Moreover, in the light of definitions just made, 
we can further identify the Markov triples, x, y and z with the regions X, Y, Z 
and W as follows : x=<f>(X),etc. where <j>(X) = tra,etc. Using Eq.s (2. 7), (2. 8) 
one can easily obtain the following result 

z,w= h(x, y) = ^(1 ± ^1-4(^ + 1)). (4.17) 
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This results prompts us to rewrite Eq.s(2.7),(2.8) in the form 



^- + ^ + ^ = 1 (4.18) 

yz xz xy 

and 

■i-+™=l. (4.19) 

xy xy 

Written in such form, these equations acquire new geometrical meaning. Given 
a directed edge e e E(J2) > define ^(e) = ^ so that for each edge e we have, 
instead of Eq.(4.19), 

ip{e) + ip(-e) = 1 (4.20a) 

(edge relation), while instead of Eq.(4.18), we have 

^(ei) + V(e 2 )+^(e 3 ) = 1. (4.20b) 

(vertex relation). It could be rather easily shown that h(x, y) > h(x)+h(y). Given 
this fact, we observe that when 2< | cc | < oo, h(x) is real. Accordingly, we ex- 
pect that |ie|, \y\ , \z\ > 2. Consider now some vertex v* G V(Y]) and consider 
the set T n (v*) which is a tree spanned by all vertices which are at the distance 
at most n from the "seed". Let Cl n (y*) be the set of all complementary regions 
meeting T n (v*) while C n (v*) be the corresponding subset of edges. Then, it can 
be shown [85] that 

£ h(<f>(X))<± ]T *& = \ ( 4 - 21 ) 

Actually, the inequality above can be replaced with equality and the above 
equality is known as McShane identity [86]. Since, in spite of its importance, we 
are not going to use it, we are not going to discuss its significance. The above 
result is mentioned only because it is of interest to inquire what happens with 
the convergence of some other functions f(X). Let, for instance, /: O — ► [0, oo) 
has a lower Fibonacci bound (to be defined shortly below), then, according to 
Bowdich [35], the series 

F(f) = lf(X)]- s (4.22) 
xen 

converges for all s>2. A lower Fibonacci bound for function /(J*f)on 0„ S fl 
exist if there is some constant k>0 such that f(X) > kF e (X) for all but finitely 
many X <G 0„. The function F e (X) can be chosen as the length L of the reduced 
word W r [ 35 ], e.g. see Eq.(2.5). Evidently, L—n in view of the results just 
presented. If this is so, the question arises: how n depends on XI Stated 
alternatively: is there way to convert the summation over X (that is over fl n ) 
into summation over nl Fortunately, the last problem can be easily solved. 
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Since X, Y, Z and W are in one to one correspondence with the Farey numbers 
all these numbers are coprime to each other [35]. Surely, that they are also 
coprime to n.This means that the length L can take any value n at precisely 
2<p(n) regions of 0„ where (f>(n) has been defined by Eq.(1.5) as number of 
numbers less than n which are prime to n. The factor of 2 is easy to understand 
if one puts n=l. Hence, Eq.(4.22) can be rewritten as 

F(f) = 2 £ 0(„)„- = 2^-^ (4.23) 

n=l ^ > 

This result coincides with earlier obtained, Eq.(4.13), thus providing an inde- 
pendent support to earlier developed thermodynamic formalism. 



5 Further developments 

5.1 The circle packing and the Farey numbers 

Rademacher [48] had found alternative geomcrtic representation of the Farey 
numbers which is worth discussing now since it has physical significance which 
ultimately goes far beyond the leaky torus model. To facilitate our reader's 
understanding, we would like to remind at this point few relevant facts from the 
theory of Mobius transformations. To begin, using Eq.(2.38) we notice that the 
point oo is the fixed point of the Mobius (in our case, modular) transformation 
as long as c = 0. If this is the case and, taking into account that ad = l,we 
obtain a transformation of the type 

z' = a 2 z + ab (5.1) 

which will have another fixed point z* — ab/(l — a 2 ) ^ oo as long as a ^ 1. 
This cannot happen, however, as we had discussed before Eq.(2.69). Therefore, 
the only possibility which is left to us is a = 1 and, hence, we are left with the 
parabolic transformation 

z' = z + b. (5.2) 

This transformation fixes infinity and makes all integers equivalent (e.g. discus- 
sion after Eq.(2.58)). It is directly associated with the presence of a puncture 
(cusp) in the case of torus as can be easily proven [87]. In H 2 model realiza- 
tion of hyperbolic space consider the geodesic (or the set of geodesies) which 
pass through point at infinity. Surely, these are just semi infinite rays which 
are perpendicular to the real axis. Accordingly, the associated horocycles are 
just set of lines which are parallel to the real axis. Consider, in particular, the 
horocycle which is located at the vertical distance 1 from the real axis. The 
parabolic transformation, Eq.(5.2), is not going to change the location of such 
horocycle. But, in general, it is known, that the Mobius transformations trans- 
form lines into lines and circles into circles. Hence, for the above horocycle 
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Figure 10: Circle packing associated with the Farey numbers. 



(actually a circle located at infinity) there must be a modular transformation 
which transforms a circle at infinity to a circle located at some point on the real 
axis. In particular, e.g. see Ref.[87], the transformation z i — ► — 1/z maps the 
horizontal horocycle at unit height to the horocycle resting at the origin and 
having diameter 1. This example can be easily generalized. To this purpose it is 
convenient to rewrite the modular transformation, Eq.(2.37), as follows 

z'-- = - 2 . 1 —/ v (5-3) 

c c 2 (z + d/c) 

Let now Z' = z'—— and Z = c 2 (z+d/c) then, we obtain back the transformation 

Z' = —1/Z. Let, now z* = Z/c 2 . By this transformation a circle of radius 1/2 
becomes a circle of radius l/2c 2 touching the real axis at zero. Since c can be 
only an integer number it is clear that all radii are less or equal to 1/2. Finally, 
by using the transformation z = z * —d/c we shift the point of tangency for 
such circle to the location z = —d/c. Consider a very special case first: c=l. 
Then, our circle is going to touch the real axis at some integer point -d of real 
axis. It is clear now that for different d's we would have different circles and 
that all these circles are going to touch each other. Moreover, and this is not 
difficult to prove [98], for c different from 1 the corresponding circles are all going 
to touch each other as it is depicted in Fig. 10. Clearly, Fig. 10 represents a very 
coarse picture since the number of circles is countable infinity between every 
neighboring pair (unit interval) of integers. Given this, one may think about the 
distribution of sizes of such circles. Following Sullivan [88], we say that two real 
numbers have the same p— size if they belong to one of the intervals (p n+1 ,p n ). 
Then, we can group the circles into collections whose diameters have the same 
p— size. The number of circles of a given size s = (p n+1 ,p n ) within a unit 
interval is the number of pairs (p, q) with p < q and p relatively prime to q (so 
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that s <~ l/q 2 ) which is again the Euler <j>(q) totient function, Eq.(1.5). This 
fact can be easily understood if we recognize, that the horocycles which are 
located at positions p/q within unit interval are having sizes of order l/q 2 . 

Following Sullivan [88] (and also Ref.[89]), let us consider definitions of pack- 
ing and covering Hausdorff measures. Adopted to our case, we have to consider 
a (closed) subset A of R 2 and to cover this subset by discs of radii n, r 2 , 
all less than some e > 0. Consider now the sum S c = Y)f V'( r i) where V'( r i) — r t 
with S being some "critical exponent" (Hausdorff dimension) of the set A. The 
covering Hausdorff ip— measure of A is the limit (as e — > 0) of the inftmum of 
S c the exponent S being a fractal dimension of A. Analogously, one can define 
the packing Hausdorff measure. In this case one should consider an open set A' 
and to cover it by the set of disjoint (that is non touching) circles. Then, one 
has to consider the supremum for the analogous sum S p . In general the packing 
and the covering Hausdorff dimensions are not the same. Such technicalities 
are needed for description of limit sets of Fuchsian (Mobius in general) groups 
with cusps (for a quick introduction to this subject, please consult our earlier 
work, Ref. [25]). Depending upon the rank of the cusp one should use either S c 
or S p [90]. At this point it is just sufficient to take into account that in our case, 
in view of the results just obtained, 

J2 (rize)*=f>(n)n-« = £^. (5.4) 

circles n—1 

Hence, we have reobtained again the result, Eq.(4.13), which now acquires com- 
pletely new meaning. It should be apparent at this point, that one still can do 
much better if one recalls all relevant facts about the Patterson-Sullivan mea- 
sure of the limit set A which we had discussed in our earlier work, Ref. [25] , in 
connection with AdS/CFT correspondence. 

5.2 The Eisenstein series and the S-matrix 

If p(x,y) is the hyperbolic distance between points x and y € H 2 then, the 
Patterson -Sullivan measure can be constructed with help of the Poincare' series 
gs{x, y) defined as 

for some Fuchsian(or Mobius, in general) group T. The factor 5 is responsible 
for convergence/divergence of gs and its threshold value is associated with the 
fractal dimension of the limit set A (which is closed (sub)set of the boundary 
at infinity, in our case, S^). According to the theorem of Beardon and Maskit 
(Theorem 5.1. of Ref. [25]) the limit set A of the discrete group T is made of 
parabolic limit points (i.e. those which are associated with cusps) and conical 
limit points (i.e. those which are associated with the fixed points of hyperbolic 
elements of T) and the conical limit points always lie outside of the cusps. 
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Figure 11: Entering and leaving the "black hole" (the cusp) while moving along 
the hyperrbolic geodesic 

This means that motion along the hyperbolic geodesies (e.g. see sections 2-4) 
is accompanied with countable infinity of events associated with such geodesic 
entering and leaving the corresponding horocycle associated with the cusp as it 
is schematically depicted in Fig. 1 1 . This fact was noticed in the paper of Sullivan 
[89] who also had estimated the characteristic time of this process (that is for 
the motion with unit speed along the geodesic (world line) one can estimate 
how long such motion spends inside the cusp). For additional illustration and 
more details, please, see Fig. 5 of our earlier work, Rcf. [25], and the comments 
associated with it. Let our geodesic begins at the point z* and ends at the point 
w* of real axis. These points are chosen in such a way that it passes through 
the point z which is the top of horocycle located at the point -d/c of real axis 
and through the point w of the horocycle located at infinity. Hence, the point 
w = i + x, the point z = iy — d/c (and y = 1/c 2 ). It is always possible to find 
a transformation g(z) such that g(z*) = 0, g(w*) = oo, g(z) = iy and g(w) = i 
[91]. In view of this, we obtain, 

p(z,w)=ln(±). (5.6) 
Using this result in Eq.(5.5) we obtain, 

£M) = 5>y)*- (5.7) 

This is just the Eisenstein series. Evidently, the subset T must correspond to 
subset of closed nonperipheral curves (that is, it does not contain an accidental 
parabolic elements) on the torus which belong to f2 (defined after Eq.(2.5)). 
Fortunately, the properties of the Eisenstein series, Eq.(5.7), are well known. 
This fact allows us not only to provide different interpretation to our main 
result, Eq.(4.13), but allows, in principle, to generalize the obtained results to 
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the Ricmann surfaces of higher genus. To understand why this is so, we would 
like to reproduce some results from our previous work, Rcf[25], at this time. 

In the upper half space model realization of d+1 dimensional hyperbolic 
space H d+1 the hyperbolic Laplacian acts on some function /(x, z) according 
to the following prescription 



A h /(x, z) = z 2 



z oz 



z > 0, x e R d . (5.8) 



In 2 dimensions the second term vanishes. We would like to keep it, nevertheless, 
since the obtained results can be immediately generalized (see below). For any 
d > 1 the eigenvalue equation for the hyperbolic Laplacian reads: 

S S 

*hz* = -(--d)s*. (5.9) 

In the case of two dimensions, of course, we have to replace z by y. If 

A h f(x) = F(x) (5.10) 

and x={x,z} then, for any 7 <G T where T is the group of isometries which leave 
jjd+i invariant, we obtain 

A h f( 1X ) = F(jx). (5.11) 

5 1 
This means that not only z ? is an eigenvalue of Ah but (72) 2 as well. Of course, 

the linear combination is also an eigenvalue of A^. Thus, not only for d=l but 

for any d> 1 we obtain, 

A h E(z,S)= 6 -(^-d)E(z,S). (5.12) 
Surely, for d = 1 we relabel z as y. In this case, following [92] we obtain, 

E(y, &)=yi+ M^y 1 ^ + 0( e - 2 ^), (5.13) 

which is sufficient for large y's. Here the scattering S-matrix f r (S) [12] is given 
by 

This result should be compared with Eq.s(4.13) and (5.4). In order to generalize 
these results we would like to mention several useful properties of the S-matrix 
</2 r (<5). Fist, it can be shown [92], that 

E{y,S) = i Pr {S)E(y,l-S). (5.15) 
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Using this result, we introduce new variable 5 — ^ = £• Eq.(5.15), when written 
in terms of this new variable, acquires the following more symmetric form 

E{y^+ l -) = ^ T {i + \)E{y^-\). (5.16) 

Let now £ — ► — £ in Eq.(5.16). By combining thus obtained equation with 
Eq.(5.16) we obtain very important relation (unitarity condition) [93]. 

Being armed with this results we can now proceed with generalizations. 



5.3 Hyperbolic Dehn surgery once again 

In this work we only provide an outline of the relevant results leaving more 
detailed discussion for future publications. To begin, let us notice that so far 
we were able to obtain all results of this paper only because the Teichmuller 
space of the punctured torus happen to coincide with the hyperbolic upper plane 
Poincare' model H 2 so that the motion in the Teichmuller space coincides with 
the motion in the hyperbolic space H 2 . This is not the case for the Riemann 
surfaces of higher genus [41] and, hence, at first sight, the results of this paper are 
not extendable to surfaces of higher genus. Very fortunately, this is not the case 
as we would like now to argue. Let us recall, that we have started our discussion 
with the figure eight and the trefoil knots in section 3. Then, we had considered 
the corresponding 3-manifolds created by the mapping torus construction and, 
after this, we had considered the incompressible surfaces (in section 3.3). Three- 
manifolds associated with these surfaces had been obtained with help of the 
Dehn surgery (Dehn filling) from the "parent" 3-manifold (e.g. see Remark 3.7) 
which is just that for the figure eight complement in S 3 . Neumann and Zagier 
in the benchmark paper [94] had developed a sort of perturbative calculations 
which allow to estimate volumes of 3-manifolds obtained from the parent 3- 
manifold by the operation of Dehn filling. In particular, for the decendants of 
the figure eight hyperbolic manifolds they obtained the following estimate for 
the volume of M(p,q) 3- manifold 



2V3tt 2 4V3(p 4 - 72pV + 144q 4 )7r 4 
p 2 + 12q 2 + 3(p 2 + I2q 2 

(5.18) 



Vol { M( P , q)) = Vol(M 8) - + — - - — + ... 



where Vol(M s ) is known to be 2.0298832... The obtained result is the simplest 
in the chain of results obtained in this reference. In general, Vol(M(p,q)) < 
Vol(M) where M is parental 3-manifold. This constitutes the essence of Thurston's 
Dehn surgery theorem [19,76]. 

In our previous work, Ref. [25], we had considered 3 manifolds with cusps, in 
particular, the 3-manifold for the figure eight knot contains just one Z©Z cusp. 
It can be shown, e.g. see the Appendix A. 3, that all hyperbolic 3-manifolds 
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associated with knots and links are Z©Z cusped: one cusp for each embedded 
circle S 1 . Hence, such manifolds are necessarily noncompact but, nevertheless, 
of finite volume. Surely, there are 3-manifolds without cusps too and these are 
related to those with cusps. According to Thurston (e.g. see section 5.33 of [19] 
), all 3-manifolds M(p 7 q) obtained by Dehn filling are without cusps. And, 
moreover, 

lim Vol(M(p,q)) = Vol(M). (5.19) 

(p,g)->oo 

Obtained result is naturally extendable to the case of k cusped manifolds. By 
introducing notations 

M k = M (pi gi ... , pk g k) 

the following general result was obtained by Neumann and Zagier [94] 

fe 

Vol(M k ) = Vol(M) - £ V Li + 0(L}) (5.20) 
2 i=i 

where Lj is the length of short geodesic r y i on which is just the length 
of the core curve of the solid torus added at i-th cusp. The result of major 
importance for us is the observation by Neumann and Zagier that the difference 
of volumes depends to a high order only on the geometry of cusps and 
not on the rest of M [94]. Given this observation, the following question can 
be asked immediately : is it possible to reobtain the results Eq.s(5.13)-(5.15) 
for cusped 3- manifolds ? The answer is YES! if these manifolds belong to the 
arithmetic hyperbolic 3-manifolds. We provide a condensed summary of results 
related to such manifolds in the Appendices A. 2 and A. 3 while in the main text 
we discuss the scattering theory for such manifolds. 

5.4 Scattering theory for arithmetic 3-manifolds 

Let is begin with observation that Eq.(5.12) is valid for any d> 1. Accordingly, 
the asymptotic expansion, Eq.(5.13), also survives and acquires the following 
form (for H 3 ) 

E(y, S)=y* + $ r (i)!/ 2 ^ + 0{e~ cv ) (5.21) 

with c being some constant. The question immediately arises: will the S-matrix 
$r(<5) have the same form as in Eq.(5.14)? In general, the answer is "NO" but 
for the arithmetic 3-manifolds the answer is "YES" as had been demonstrated 
by Sarnak[26]. It is instructive to restore some of his calculations now. In order 
to do so, we have to think about the ways the group PSL(2,Z) can be extended 
in order to be used for description of 3-manifolds associated with the punctured 
torus bundle. 

To this purpose we notice that the group PSL(2,Z) is just a subgroup of 
PSL(2,R) which is group of isometries of H 2 .The group of isometries of H 3 is 
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PSL(2,C) as is well known [25]. Hence, we have to figure out what kind of 
subgroup of PSL(2,C) is analogous to PSL(2,Z). This problem was actually 
solved for the torus bundles in Ref. [95] and for surfaces of higher genus in very 
important paper by Margulis [30] to be discussed further in Appendix A. 3. For 
the torus bundles, naturally, one should try to find solutions of Eq.(2.7) in the 
complex domain. In Ref. [95] it was shown that the complex counterpart of the 
triple (3,3,3) is 



and the rest of triples also can be obtained . The numbers above are integers 
which belong to the ring of integers O3 in the quadratic number field Q(\/— 3) 
,e.g. read Appendix A. 2. Hence, for the punctured torus bundle the analog 
of PSL(2,Z) is PSL(2,03) and the punctured torus fiber bundle is 3-orbifold 
H 3 /PSL(2, O3) as was proven rigorously in Ref. [95]. The question remains: 
can this result be generalized to the noncompact (that is having punctures 
(or cusps)) 3-manifolds originating from automorphisms of Riemann surfaces of 
genus higher than one? The answer is "YES", provided that discrete group of 
isometries is Bianchi group PSL(2,0d), where 

-d = 1, 2, 3, 5, 6, 7, 10, 11, 14, 15, 19, 23, 31, 35, 39, 47, 55, 71, 95, 119. (5.23) 

This remarkable result can be found in Ref. [96] . What is even more remarkable, 
that it is just a refinment of the very comprehensive work by Bianchi [28] which 
was completed already in 1892 ! In Ref. [97], Proposition 1, Reid had proved the 
following theorem 

Theorem 5.1. Every non- compact arithmetic Kleinian group is conjugate in 
PSL(2,C) to a group commensurable with some Bianchi group PSL(2, Od)- 

Remark 5. 2. a) The definition of arithmcticity is rather involved and is dis- 
cussed in the Appendices A. 2 and A. 3. The above theorem is complementary 
to the theorem by Margulis ( e.g. Theorem A. 3. 11 of Appendix A. 3.) which 
is much more comprehensive. b) The noncompactness implies existence of the 
parabolic peripheral subgroups (e.g. see section 2).c)The notion of commensu- 
rability between groups can be found, for example, in Ref. [98], and is associated 
with set theoretical definition of overlap (i.e. n) between sets. 

Corollary 5.3. Every cusped arithmetic hyperbolic 3-orbifold is commensu- 
rable with Bianchi orbifold H 3 / PSL(2, Od) where Od is ring of integers in the 
quadratic number field Q( \J — d). 

Remark 5.4. a)There is one-to-one correspondence between the rings S 1 in 
somc(arithmctic) link and the cusps (e.g. read AppendixA.3. for more details), 
b) The non-compactness is caused by cusps. c) Cusped 3-manifolds do have finite 
volume(e.g. read section 5.3 again), d) The arguments of section 5.1. strongly 




(5.22) 
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suggest that 3-manifold associated with the figure eight knot is arithmetic. This 
is rigorously proven by Reid in [97] who also proved that the figure eight is the 
only knot which is arithmetic.e) Corollary 5.3 and Eq.(5.23) suggest that there 
is a countable infinity of arithmetic links. 

Being armed with these results we can now discuss the results of Sarnak 
[26]. To this purpose, we have to rewrite Eq.(5.7) in a more convenient (for the 
present purposes) form. We have (z = x + iy): 

oo oo 

E(y,6) = yUJ2 E E 7 71 ( 5 ' 25 ) 

c=l (c,d)=l m=-oo U c ( x + m ) + dof + C 2 y 2 

e.g. see Lax and Phillips [12], pages. 171-172. Here in the second sum summation 
is over d provided that c and d are relative primes, i.e. (c, d) = 1, d=do+mc 
and 0< g? <c. The summation procedure is explained in great detail in the 
same reference. In the Sarnak's case one needs only to replace x by z and y 
by t (t >0) (in the H 3 model realization of hyperbolic space). The details of 
calculations can be found in Rcf.[26] so that the final result, indeed, has the 
form given by Eq.(5.21) where for the quadratic number field Q(\/— d) of class 
number one(e.g. read Appendix A. 2. for explanation of terminology) $r(<5) is 
given by 



,,,, ( ,j : - 11 { r> 2(i) 

with V(Fl) being defined as 



provided that 



, s/D/2 if D = 3(mod4) 
V ^~< V^ifD^3(mod4) (5 ' 27) 



d = 



-D ifD=3(mod4) 
AD if D ^ 3 (mod 4) 



Since, according to Appendix A. 2, Eq.(A.2.12), the Dcdckind zeta function ( d (s) 
has the same pole as the ordinary zeta function, the result, Eq.(5.21), can be 
used instead of earlier obtained Eq.(5.14). This time, however, one can get much 
more. In particular, one can get the volume V(Fjj) of the associated cusped 
hyperbolic 3-manifold. The expression for the volume is known to be [19,26,98] 



V{F D ) = W A ^ {2) - (5.28) 

This result can be obtained directly from the Eisenstein series, Eq. (5.25). It 
comes as the residue at the pole | = 2 of the Dcdekind zeta function in the 
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numerator Eq.(5.26) [26]. That is 



Res[E(y,6), 6 -=2} = ^±. (5.29) 

Since both the residue and V(Fl) are known, the volume, Eq.(5.28), is obtained 
using Eq.(5.29). 

Consider now generalization of the obtained results to the case of multiple 
cusps (that is to the arithmetic links (e.g. see Appendix A. 3). Following Efrat 
and Sarnak [27 ],and also [98], consider a complete set { K\ = oo, k 2 , K n } of 
the nonequivalent cusps and let Ti be a subgroup of T C PSL(2, C) for which 
Ki is the fixed point that is I\ is stabilizer in Y of the i-th cusp. It is always 
possible to select transformation /^(/t,) = oo where p i £ T.Then, 

Pi V iP T'=^ \ ^ (5.30) 

which coincides with the matrix U\ defined in Eq.(A.3.1) with uu being a unit of 
quadratic number field. For each cusp one can define its own Eisenstein series 
analogous to Eq.(5.7). Also, by analogy with Eq.(5.21) we write now 

Ei(y,8) =5 ijy t +$ i , j (% 2 "5 , i,j=l,...,n . (5.31) 

It can be shown, that the S— matrix obeys the matrix equation 

E(y,6) = $(S)E(y,2-S) (5.32) 

analogous to Eq.(5.15) discussed earlier. Finally, the analog of the Eq.(5.14) is 
the determinant of the matrix <!>((>) given by 



Ah( 5 



det<f>(6) =cons V ny2 ' (5.33) 



to be compared with Eq.(5.26). £,h( s ) f° r the quadratic number field of class 
number one is given by 

^( S ) = (Vd/2^) S r( S )C H ( S ) (5.34) 
with C H (s) being defined by 

(h(s) = 1[L(s,x). (5.35) 

x 

It should be noted, that the above product over characters (which are specific 
for each quadratic field) always starts with \ = 1 • This produces L(s, 1) = £(s) 
where, again, C(s) is usual Ricmann zeta function. Hence, again, the singu- 
larities of the determinant are just those for the ratio of zeta functions as in 
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Eq.(5.4). Hence, the exact partition function for 2+1 gravity is given by the 
determinant det$, Eq.(5.33), of the scattering S-matrix 
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Appendix A.l 



Phase transitions in 2+1 gravity: analogy with Bose condensation 

Although in section 4 the limiting value of the partition function, Eq.(4.13), 
is obtained in a systematic way [9,10], the extraction of useful information from 
this result is plagued by some serious problems as was noticed already in earlier 
work of Cvitanovic [15 ]. The result, Eq.(4.13), is obtained as limiting value of 
the Farey spin chain partition function whose size tends to infinity. Since in Sec- 
tion 5 we generalize the punctured torus results to surfaces of higher genus, the 
analogy with spin chain cannot be straightforwardly extended. Hence, we need 
to develop methods of extraction of useful information without being dependent 
on spin chain analogy. 

The main difficulty with Eq.(4.13) lies in the fact that it contains singulari- 
ties. This can be easily demonstrated if we notice that for large values of n the 
Euler totient function <j>(ri) <~ n [45]. This observation allows us to write 

Z{fi)<W-\). (A.1.1) 

Since ((/3 — 1) is singular at [3 cr = 2 one is faced with the problem of removing 
this singularity in physically acceptable way. There are many ways of doing so 
but it is not our purpose here to provide a complete list of possibilities. Rather, 
we would like to notice the following. First, by direct numerical calculation one 
observes that for large (3' s the partition function Z{(3) approaches one. So that, 
indeed, for low "temperatures" the free energy is zero in accord with earlier 
results [9,10]. The authors of these papers claim that Z{(3) is zero as soon as 
(3 > 2. This can be understood only if one considers finite number terms in the 
corresponding zeta functions and takes the thermodynamic limit using Eq.(4.8). 
To be consistent, we need to do the same for (3 < 2. This is not a simple task 
however since, according to Ref.[10] (e.g. see p. 428), "For high temperatures 
(f3 < 2) the analytic continuation of Z(f3) cannot be directly interpreted as the 
partition function of the infinite chain" . This statement is rather pessimistic, 
especially, if one is thinking about extension of obtained results to surfaces of 
higher genus. Therefore, we offer here a somewhat different interpretation of 
the obtained results. 

To this purpose, let us introduce new function g^{a) defined by the following 
equation: 

oo 

9 ^ = ^=m h^p^i- (A - L2) 

71-1 

Such type of functions are well known from the theory of Bosc-Einstcin (B-E) 
condensation [99] .The analogy, unfortunately, is not complete since in the case of 
B-E condensation the exponent (3 — 1 is never zero or less than zero. Evidently, 
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for a = 1 we obtain back the Riemann's zcta function. The parameter a is 
related to the chemical potential /i in a usual way. For /i = 0we have a = 1. This 
value of a is associated with singularity which has physical interpretation. In 
properly chosen system of units equation which connects the chemical potential 
with the particle density p — y in B-E case reads [99]: 

P = 9 l {a) + VT^ (A ' L3) 

where V is volume of the system and N is the total number of particles. The 
second term in the r.h.s. represents the density of Bose condensate associated 
with zero translational mode. For finite volume this density is infinite for a = 1. 
Therefore, as usual, it is assumed that both the volume and the density tend 
to infinity in such a way that the ratio is finite number. In our case, we do 
not have the luxury of having volume term but we can extract the singularity 
in a manner similar to the B-E gas. Taking into account Eq.s (A. 1.1) - (A. 1.3) 
we can subdivide the domain of integration into two parts : a) 0< x < 1 and 
b)l< x < oo. In the first case we obtain, 

l 

h=f dxx?-* = j-^ (A.1.4) 
o 

which is expected singularity of the zeta function at already known critical value 
[3 cr = 2. As for the second part we obtain (close to criticality) 

h = j dx^ * J dx-±- x +W-2)J dx^ + ... (A.1.5) 
l l l 

Combining Eqs.(A.1.4) and (A.1.5) we obtain as well, 

/ = i(l + c l£ )(l + c 2 (/3)e 2 + ...), (A.1.6) 

where e = (3 — 2 and Ci,C2(/3) are known constants. In particular, 0.373 < c 2 < 
0.5438 for 1< f3 < 2. Taking into account that T(f3) is nonsingular in the range 
of 0s which is physically interesting, we can subtract the singular part along 
with well behaving regular part in order to get 

-W) =ln(l + c 2 (l3)e 2 + ...) > c 2 (/3)£ 2 . (A.1.7) 

This result coincides almost exactly with the estimate, Eq.(25), obtained in 
Ref.[9] (the constant c 2 is w 0.25 in Ref.[9]) for the free energy in the high tem- 
perature "disordered" (or pseudo-Anosov) phase. Hence, the rest of arguments 
of Ref.[9] now can be used. In particular, we have 

2 

-/W) - f U(x)dx < U((3){2 - [3). (A.1.8) 
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Using this result, we obtain, 

U(J3)>-- W)/(2-/3)>c 2 (/3)(2-/3), 1 < /3 < 2, (A.1.9) 
in complete accord with Eq.(4.15) of the main text. 

Remark A 1.1. a) The analogy with Bose gas condensation can be strength- 
ened, perhaps, even more if one notices that the criticality condition (3 cr = 2 
coincides with the simple estimate made in famous Kostcrlitz and Thouless pa- 
per [23]. In this paper the Coulomb gas model is used for description of phase 
transitions in two dimensional liquid crystals, liquid helium, etc. This, in part, 
explains the success of such type of models prevailing so far in physics literature 
[20]. b) In the case of liquid helium Feynman [100] had obtained equation very 
similar to Eq.(A.1.3) and Landau and Lifshitz , e.g. see section 27 of Ref.[101], 
also resort to a simple subtraction of the undesirable singularity at zero tem- 
peratures. Hence, intuitive physical arguments may provide a reliable guidance 
needed for further refinements of the results just presented. 

Appendix A. 2 

Some results from the algebraic number theory 

Theory of arithmetic hyperbolic manifols has strong connections with the 
algebraic number theory [102,103]. In this appendix we provide a condensed 
summary of ideas on which such theory is based. Selection of topics is, naturally, 
subjective and is meant only to encourage interested reader to read much more 
comprehensive texts. 

Development of number theory is associated with the desire to extend a 
simple idea, known to nonprofessionals, that every nonnegative integer n can 
be uniquely represented as 

n=p?---p? (A.2.1) 

where pi is some prime number and on > is some integer. This result becomes 
much less obvious if one is willing to extend it, say, into domain of complex or 
algebraic numbers. For instance, 

6 - 2 • 3 = (1 + V^5)(l - = (4 + Vl0)(4 - \/To). (A.2.2) 

Such non uniqueness is highly undesirable since many theorems of usual arith- 
metic brake down. Solution of the nonuniqucness problem is one of the major 
tasks of the modern number theory. 

Let us define a rational number £ =p/q as a root of the equation 

q£-p = (A.2.3) 
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with p and q being some integers. An integer is then a root of an equation 
with first coefficient q=l. Analogously, an algebraic number £ is any root of 
an algebraic equation 

a n x n + a n ^ix n ^ 1 H h a = (A.2.4) 

where the coefficients a n , ■ ■ -ao are rational integers. Accordingly, an algebraic 
integer is any root of the monic polynomial, e.g. see Eq.(3.22). 

Remark A 2.1. Since only fibered knots and links are relevant for dynamics 
of 2+1 gravity, e.g. read section 3, and since such knots/links are associated 
with monic Alexander polynomials, we conclude, that the arithmeticity is an 
intrinsic property of 2+1 gravity. Surely, much more is associated with this 
observation as we shall demonstrate shortly. 

In particular, let us consider the module (that is the set which is closed 
under operations of addition and subtraction) of quadratic integers, i.e. those, 
originating as roots of quadratic equation with 02 = 1. We had discussed the 
roots of such equations already, e.g. see Eq.s(3.2), (3.9), (3.18)and (3.35). Now 
we know that all these roots are quadratic integers. As with usual (rational) 
integers, it is important to define the units, especially if we would like to go 
from the module to a ring ( which is module where the multiplication operation 
is defined which keeps all the members of the set within the set) or to a field 
(where, in addition, the operation of division is defined). It is clear that the 
typical representative of a quadratic integer £ should look like 

e= a + bV ^ (A.2.5) 

c 

where a, b, c and d are some integers. Moreover, more detailed analysis shows 
that d can be chosen as square free. Let us notice, that in the case if d=-l, we 
obtain more familiar field of complex numbers. Geometrically, the integer lattice 
in such field can be built using the basis [1, i]. Therefore, in general case we can 
think of a basis [1, u] where 

u=^(d + Vd) (A.2.6) 

with d=B if D=2,3 or rf=4D if D=l (mod 4) e.g. see Eq.(A.2.8) for explanation 
of notations. The square free number d is called discriminant of the field. This 
result confirms, in particular, that the Markov triple in Eq.(5.22) is indeed made 
of integers since |(3 — \/— 3) is just the conjugate of ^(3 + \/— 3). For different 
d's one will have different rings of quadratic integers. To make rings, once the 
integers are defined, one should also define the units. This can be done again by 
analogy with complex analysis. Using Eq.(A.2.5) we define a norm N(£) as 

N(Q=t-t (A.2.7) 
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where £ is conjugate of £, i.e. the same as £ in Eq.(A.2.5) but with b — > —6. In 
the number theory the norm can be both positive and negative. For quadratic 
imaginary field it is positive, however. By definition, units are such numbers 
whose norm is unity. In the case of usual complex numbers, one has therefore 4 
units: ±1 and ±i , since in both cases the norm is one. In the case of d~-3 ring 
the units are ±1, ±p, ±(1 + p) and p 2 where p = ^(— 1 + \/^3),in all other cases 
the number of units is 2, e.g. see pagel56 of Ref. [102]. The role of units can be 
easily understood by analogy with more familiar case of ordinary integers . That 
is if a and b are integers then, a = 1 • b = b ■ 1. Units play the same role for the 
ring of quadratic numbers. These numbers are called associates if they differ by 
unit factor. Next, one defines primes as those quadratic numbers whose norm 
is a rational prime. In particular, we notice, that the Markov triple, Eq.(5.22), 
is made of primes. Once the primes are defined, one is naturally interested in 
finding an analogue of the Eq.(A.2.1) for quadratic integers. This happens to be 
not a simple task: not all d's will lead to unique decomposition into primes. For 
reasons which will become clear upon reading of Appendix A. 3, in the case of 
gravity only negative d's are relevant. Among those which allow unique prime 
decomposition are: -d=l, 2,3, 7, 11, 19,43, 67 and 163 [102,103] to be compared 
with Eq.(5.23). Such comparison indicates, that some d's listed in Eq.(5.23) 
suffer from non uniqueness. Such non uniqueness, had lead number theoreticians 
to introduce the concepts of ideals and class numbers. Class numbers reflect the 
extent to which integers of the field deviate from uniqueness of decomposition. 
Naturally, if the decomposition is unique, the class number is one. The ideals 
can be defined set theoretically and geometrically as follows. Suppose,we have 
some set I and a larger set L, then V a, [3 e I and V£ e L we have 

a + /3 e I (module property) 

and 

a£ G / (ideal property). 

Geometrically, this can be understood as follows. Consider, without loss of 
generality, two dimensional lattice £ built on vectors [l,u;] then, for any rational 
integer n the point [n, nu>] belongs to £. Now one can think of different classes 
of rational integers and, hence, of different sublattices 9Jt C £ . This is easily 
accomplished by introducing the residue classes modulo m. Recall, that notation 
(equivalence relation) 

x = y (modm) (A. 2. 8) 

means that x-y=mk where k is some integer. The above notation is just the 
statement that when both x and y are divided by m they both have the same 
residues. In particular, the statement x= O(modm) means that m divides x. 
Let m have the same presentation as n in Eq.(A.2.1), then any number x rela- 
tively prime to m may be determined (modm) by the equations x= x^modp^*) 
, (xi,Pi) = 1, i=l, 2,..., k. The number of residue classes is just the Euler 
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function </>(m). Now, if for a given sublattice DJl and vectors x and y e £ 
we have equivalence relation 



x ee y (modM) 



(A.2.9) 



this means only that x-y€ 9Jt. The number of different sublattices within given 
lattice (i.e. the number of different residue classes for ideals) is called index 
j. Evidently, j=[£/9Jt]. Equivalently, it is also called a norm of an ideal £ 
so that j=N[£]. ft happens, that in the case of quadratic fields calculation 
of the norm is relatively easy [103]. If a is quadratic or rational integer, then 
N[£]=a 2 =|A r (a)|. For both ideals and integers the respective norms obey the 
following composition law 



where the Gothic letters stand for the ideals and the italics for numbers. The 
question arises: if there are prime numbers are there prime ideals? The answer 
is "Yes" and since prime numbers play the central role in theory of rational 
integers, one should expect that the prime ideals play the same role in the 
number theory This is indeed the case, but, unlike numbers, every ideal is 
decomposable into prime ideals uniquely. This is the main reason why these 
objects were introduced. 

Remark A 2.2. In the case of "quadratic fields" there is one to one cor- 
respondence between the "quadratic ideals" and the quadratic forms[102 ]. We 
deliberately omit discussion of this fact to save the space referring interested 
readers to literature [102-104]. 

Without going into details of such decomposition, we would like to take 
advantage of the composition laws, Eqs. A. 2. 10-11. To this purpose we need the 
following 

Theorem A 2.3. The rational prime p factors in the quadratic number 
field Q(-\/d) (irrespective to the sign of d), with accuracy up to 
multiplication by unit, as follows: 

a) p does not factor (p is inert in the field d) if and only if (d/p)=-l, 

b) p splits in d into two different factors (e.g. z and z' which belong to 
QfVd)) if and only if (d/p)=-l, 

c) p is ramified in d (i.e. p~z 2 ) if and only if (d/p)=0. 

In the above cases (d/p) is the Kronecker symbol which sometimes is denoted 
as x(p) for fixed d. x(p) is actually a character (the Dirichlct's character to be 
exact) of an Abelian group, e.g. of £/9JT . This symbol should not be confused 



N(o)N(b)=N(ob) 



(A.2.10) 



and 



N(a)N(b)=N(ab) 



(A.2.11) 
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with Sij known in physics. At the same time, the above theorem may serve as 
a definition of such a symbol. We shall adopt this point of view in this work. 
Then, by analogy with the theorem just stated, one can formulate the analogous 
theorem for the ideals [103] 

Theorem A 2.4. The quadratic prime ideals (p) are related to integers p 
in the rational field in the following possible ways: 

a) (p)=(p) , or (p) does not factor, i.e. inert, then N (p)=p 2 ; 

b) fp)=pip 2 , or (p) splits, then Nfpi)=Nfp 2 ) = p; 

c) (p)-pl, or (p) ramifies, then NYpi) =p 

Being armed with such results, we are ready to introduce the zeta function 
for the field whose discriminant is d: 

Us) = E [ N ( fl )P = - MpTT 1 (A.2.11) 

a p 

where the last product is over all prime ideals. In view of the Theorem A 2.4., 
this result can be rewritten in more familiar terms as follows: 

C» = CM £ ^ = C(s)L(s; d) (A.2.12) 

n— 1 

where L(s,d) is known in the literature as the Dirichlet L-series. It can be shown 
[102-104], that such defined zeta function is having a pole singularity only at 
s=l, i.e. the singularity only comes from the usual zeta function. The residue 
of the pole at s=l is much more interesting in the present case as it is explained 
in the main text in section 5. Finally, we need an ideal analogue of the Euler 
totient function <f>(n). To this purpose we take into account that an analogue of 
Eq.(A.2.1) is the same equation in which all numbers are replaced by the ideals. 
Then, for an ideal a we have 

*(«) = N(a)n(l-i). (A-2.13) 

p/a 

where p runs through the distinct prime divisors of o . 
Appendix A. 3 

Connections between the knot theory, theory of hyperbolic 
spaces and the algebraic number theory: implications for 
2+1 and 3+1 gravity 

For the sake of saving space and to avoid repetitions, we expect our readers 
to be familiar with our earlier work, Ref.[25], while reading this Appendix. Let 
us begin with the group PSL(2,C) of isometries of hyperbolic space H 3 . Only 
discrete subgroups of PSL(2,C) known as Klcinian groups are of interest [98]. 
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The group PSL(2,C) is just the projectivized version of SL(2,C). The following 
theorem for SL(2,C) is of particular interest, e.g. see Ref.[98], chr.l,. 



Theorem A. 3.1. The group SL(2,C) is generated by two elements 

Ui=(q I ) and U 2 = ( J "q 1 ) , where a e C. (A.3.1) 



Remark A. 3. 2. This theorem is an analogue of the well known fact [105] 
that the group PSL(2,Z) can be generated by the same two elements (upon 
projectivization) with a=l. Hence, the parabolic elements play a very special 
role in PSL(2,C) too. 

The parabolic subgroups of PSL(2,C) are abelian subgroups associated with 
Z and Z0Z cusps. If T ePSL(2,C) contains an abelian subgroup of rank one, 
i.e. Z-cusp, then the associated fundamental domain for T is non compact and 
of infinite volume [29]. Hence, of no interest to us (additional reasons are given 
below). If, however, T ePSL(2,C) contains an abelian subgroup of rank two, 
i.e. Z©Z cusp, then the group T is not cocompact but of finite covolume. The 
following theorem was used in our earlier work, Ref . [25 ] , and can be stated now 
as follows : 

Theorem A. 3. 3. Let T ePSL(2,C) be a discrete group of finite covolume, 
then r has only finitely many classes of cusps. 

The proof can be found in Ref. [98], please, read Chr.2, Proposition 3.8. .To 
insure that the subgroup T is discrete the following theorem is the most helpful 
[98] 

Theorem A. 3. 4. (Jorgensen's inequality) If <A,B>is non- elementary dis- 
crete subgroup of PSL(2,C), then, the following inequality holds 

\tr 2 A-A\ + \tr[A,B] - 2| > 1. (A.3.2) 



Remark A. 3. 5. a) In the case of punctured torus extension of Markov 
triples Eq.s(2.7),(2.8) to the field of quadratic integers was discussed in Ap- 
pendix A. 2 and in the original works, Rcfs. [35,95]. Such an extension is in 
accord with Jorgensen's inequality b) This observation provides us with an 
easy way of definig the arithmetic Klcnian groups and, hence, the arithmetic 
3-manifolds. To this purpose, we need still one more theorem [95]. 

Theorem A. 3. 6. Let A 1 ,A 2 ,...,A n e SL(2,C) be such that A x is not 
parabolic and <Ai,A% > is not elementary. Then, the A\, A 2 , A n are deter- 
mined up to conjugacy in SL(2,C) by Tr A i: TrAiAj and TrAiAjAk with i,j,k 
e{l,2,...,n}. 
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Remark A. 3. 7. The arithmetic Klcinian groups can be regarded as im- 
ages in PSL(2,C) of subgroups of SL(2,C) so that the trace of an element in 
the Kleinian group is defined up to a sign. The traces Tr Ai belong to the 
field of quadratic integers with negative discriminant d. The negativity of d is 
explained shortly below. 

Following [106], page 192, Theorem E. 4. 4., and also Riley [29] consider com- 
plement of a link L in S 3 (which is just a union of knots C\, £fc).To this 
purpose, let Tj be an open tubular neighborhood of the knot (that is Tj is a 
solid torus with boundary removed). Then, the following theorem holds: 

Theorem A. 3. 8. Given a compact 3-manifold N there exists a link C in S 3 
which is union of knots C\, ...,£k, such that N is obtained from the manifold 

M = 5 3 \( (J T *)' 

»=l,...,k 

where T^s are pairvise disjoint open tubular neighborhoods of the knots £i,by 
the operation of Dehn filling. 

Corollary A. 3. 9. (Riley and Thurston [19,29]) // M is to be hyperbolic, 
then every non-cyclic abelian subgroup (of rank two) of iriM must be peripheral 
and correspond to the component d of link C. 

Remark A. 3. 10. Since the fundamental group iriM can be embedded 
into PSL(2,C) and since the non-cyclic peripheral (e.g. section 2 for definition) 
abelian subgroup of rank two corresponds to Z®Z cusp there is a one to one 
correspondence between the cusps and the link components. E.g. in the case 
of figure eight knot we have hyperbolic 3-manifold with just one cusp [19,29]. 
Moreover, and this is the most important, the following theorem of Margulis 
[30] requires such cusped 3-manifolds to be arithmetic. 

Theorem A. 3. 11. (G. Margulis) Let S be Riemannian symmetric noncom- 
pact space of rank greater than one and A is irreducible discrete group of motions 
of S such that S/ A has finite volume but noncompact. Then, A is the arithmetic 
subgroup in the group of motions of space S. 

Definition A. 3. 12 ( Helgason [31], Chr.5 ) Let S be the Riemannian glob- 
ally symmetric space. The rank of S is the maximal dimension of a flat, totally 
geodesic subspace M. The subspace M is geodesic at p £ M if each S-geodesic 
which is tangent to M at p is a curve in M. The submanifold M is totally 
geodesic if it is geodesic at each of its points. 

Remark A. 3. 13. To use the above definition, the following alternative 
interpretation of rank is helpful. Following Besse [32], let Q=t)+m be decompo- 
sition (into irreducible parts) of the Lie algebra of the group of motions of space 
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S where f) is non and m is abelian parts. Then, the rank is the total dimension 
of the abelian subalgebra m of g. 

Remark A. 3. 14. Since H 3 is symmetric space, Margulis theorem implies 
arithmeticity of hyperbolic 3-manifolds with Z+Z cusps. In the case of one cusp 
the complement of figure eight knot hyperbolic 3-manifold is indeed arithmetic 
as it was shown independently by Reid [97]. The arithmeticity can actually 
be anticipated based on the results of our section 3. Indeed, since surface 
automorphisms are associated with the roots of the monic Alexander polynomial 
whose coefficients are rational integers, these roots lie necessarily in the field of 
quadratic integers. The Scifcrt-fibcred phase is associated with units of the 
field while the pseudo-Anosov with nonunit integers. Since these stretch factors 
are related to the traces of the corresponding matrices, e.g. see Eq.(2.15), 
connection with Theorem A. 3. 6. is clear. 

Remark A. 3. 15. Since Margulis theorem is not limited to groups of mo- 
tions of H 3 its application can be extended to 3+1 gravity. According to 
Besse,Chr 10 of Ref.[32], every irreducible symmetric space is Einstein space. The 
classification of irreducible symmetric spaces was made by Helgason [31] and is 
also cited in Bcsse's book. From this classification it follows that most of such 
spaces have rank higher than one and, hence, should be arithmetic. It would 
be interesting to know if there is an analogue of Theorem A. 3. 3 for Einstein 
spaces. This is of interest because, from the discussion we had in section 5, it 
follows that the cusps act very much like black holes, e.g. see Fig. 11, (except in 
2+1 case the "particle" cannot be trapped by the "black hole") and, therefore, 
the role of black holes in the Universe acquires new meaning: they make our 
Universe arithmetic. If this is the case, then the Universe becomes something 
like a crystalline solid which potentially can undergo phase transitions between 
lattices of different symmetry. The symmetry of the lattice is determined by 
the discriminant of the quadratic number field. 

Remark A. 3. 16. The connection between the cusps and the black holes is 
not purely visual, e.g. Fig. 11, or formal. It is real, at least for 2+1 Euclidean 
black hole ! Indeed, following Carlip, e.g. see p. 50 of Ref[22], by the appropriate 
choice of coordinates and units the metric of the Euclidean black hole becomes 
that of H 3 , i.e. 

{ds)2 = (dxf + {dyf + {dzf ^ 

provided that in addition the following identification is made, 

(x,y,z) <~ e 27Tr+ (x cosh 2nr + — y sinh 2irr + , y cosh 2nr + — x sin nh2nr + , z) , 

which is easily recognizable as Z©Z cusp [25]. 
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Remark A. 3. 17. The reasons for existence of black holes in nature were 
discussed in the past [107 ] . Ironically, the same author had came to the opposite 
conclusions in [108], that is the black hole singularity is fictitious. The same 
conclusion is reached in Ref.[109] which claims that two dimensional quantum 
black holes are nonsingular. 

Now we have to explain why this discriminant should be negative. For 
one thing, if it would be positive the norm of the quadratic field would be 
negative. This is unphysical in view of Remark A. 3. 14. There are other reasons 
in addition which wc would like now to explain. Let us begin with famous 
example by Riley [110] who considered the mapping of the fundamental group 
7r 8 of the complement of figure eight knot into PSL(2,C).The fundamental group 
7r 8 is free non abelian group composed of two generators which has the following 
presentation 

7Tg =< X\, X2\ WX\ — X2W > (A. 3. 3) 

where w = x± X2X1X2 ■ The homomorphism of 7r 8 into PSL(2,C) can be 
achieved with help of the Theorem A. 3.1. as follows: 

x x -> A = ^ J J ^ and x 2 -» B = ( J ^ (A.3.4) 

where u> is some number, different from 1 and to be determined below. The 
matrix B is just a conjugate of U2 in Eq.(A.3.1), i.e.B = U1U2U1 with a ^ u. 
Since these matrices should obey the relation wx\ — x 2 w, replacements x\ — > A 
and X2 — > B produce the following equation for u 

w 2 +w + l = (A.3.5) 

with solution 10 = |(— 1 + \f— 3). One can easily recognize the unit p of the 
quadratic field with discriminant 3 (from Appendix A. 2). Hence, figure 8 knot 
is indeed arithmetic. From this example it should be clear, that the arithmeticity 
is associated with our choice of u. 

Consider now the subgroup of PSL(2,C) which fixes j (i.e. in the upper space 
model H 3 with coordinates z+jt (z=x+iy, t>0) we are looking for the Mobius 
transformation M(j)=j). It is easy to prove [98] that such transformation is 
associated with SL(2,C) matrix of the type 

m=(_;_ 1) (A.3.6) 

such that det M = \x\ 2 + \y\ 2 = 1 where \x\ 2 = xx and x is a complex conjugate 
of x, etc. Such matrix M is in one to one correspondence with the field Ti. of 
quaternions q known in physics literature as Hamiltonian quaternions [111] . 
Such quaternions usually are represented as q=a +bi +cj +dk where a,b,c,d are 
real numbers and k=ij. In terms of quaternions Eq. (A.3.6) can be written as 

(* + + * T-i)- ^ 
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Surely, if we write i=y/—l , then the above matrix represents the most general 
matrix associated with the unit of the imaginary quadratic field with discrim- 
inant d=-l. Naturally, one can therefore think of the most general matrix 
associated with units of imaginary quadratic fields with discriminant d other 
than -l.To this purpose one needs to extend the notion of quaternions known 
in physics. 

Let K be the number field of characteristic different from 2 (the characteristic 
p of the field K is associated with subdivision of the field into residue classes Z p , 
so p^ 2 and p is prime) and let a, b <GK be two non-zero elements. A quaternion 
algebra 7i(a,6;K) over K is generated by elements i and j satisfying i 2 = a, 
j 2 = b and ij = —ij. The elements 1, i,j, ij form a K-basis of H(a,b;K) as a 
vector space. The pair (a,b) is called the Hilbert symbol for H . In particular, 
for Hamiltonian quaternions the field is R and the Hilbert symbol is (-1,-1). If 
L is the field extension of K with y/a and y/b £ L then, repeating all the steps 
leading to Eq.(A.3.7), we obtain, 



Newman and Reid [112] proved the following major 

Theorem A. 3. 18. Let T be a non-compact Kleinian group of finite covol- 
ume. Then: 

a) r is arithmetic if and only if the trace field (e.g. see Theorem A. 3. 6) 
K = Q(^/^d) for some square free de N (e.g. see Eq.5.23) and trT consists of 
algebraic integers. 

b) T is derived from the quaternion algebra H(a,b;K) if and only if trT C 
Od for some d. 




(A.3.8) 
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